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ABSTRACT 


The  resonance  phenomenon  which  occurs  when  a  high 
•peed  turbulent  jet  impinges  on  a  large  flat  plate,  is  ex¬ 
perimentally  studied.  The  far  field  noise  is  also  inves¬ 
tigated  for  a  wide  range  of  impinging  jet  speeds.  The 
study  involves  far  field,  near  field,  and  surface  pressure 
measurements. 

A  feedback  of  pressure  waves  between  the  nozzle  and 
the  plate  is  found  to  be  the  mechanism  responsible  for  the 
resonance.  The  feedback  loop  is  formed  by  downstream 
travelling  large  coherent  structures,  and  upstream  acous¬ 
tic  waves  propagating  in  the  nearby  stagnant  region  around 
the  Jet.  The  upstream  acoustic  waves  are  generated  by  the 
impingement  of  large  scale  structures  on  the  plate.  Near 
the  norzle  exit  these  waves  are  locked  in-phase  with  the 
oscillations  of  the  thin  shear  layer,  which  is  forced  to 
roll  into  large  eddies  at  the  resonance  frequency  f  . 
Furthermore,  the  experimental  data  revealed  the  following 
interesting  results: 

(i)  The  total  period  of  the  feedback  loop  is  an  in¬ 
teger  N  times  the  resonance  period,  N  being  constant 
over  a  resonance  frequency  stage.  A  jump  from  one 
frequency  stage  t>  another  (occuring  at  certain  . 
no  e  zl  e-  to- pi  a  t  e  p«o~o,ti  on  distances),  takes  place 


through  a  quantum  unit  change  In  N,  ao  aa  to  preserve 
the  phase  lock  between  the  two  oscillations  at  the 
nozzle  . 

(11)  As  the  wavy  shear  layer  evolves  into  large  vor¬ 
tices,  a  strong  inte r action  develops  among  the  Initial 
Instability  vortices.  This  phenomenon  is  referred  to 
as  "collective  interaction",  and  is  believed  to  be  the 
mechanism  responsible  for  the  sharp  drop  in  the  pas¬ 
sage  frequency  of  the  vortices  near  the  nozzle  exit. 
The  frequency  drops  from  the  high  instability  frequen¬ 
cy  f^^  to  that  of  the  resonance  frequency  f  ^ .  Visual 
observations  in  a  forced  two-dimensional  free  shear 
layer  by  Ho  and  Huang  (1978),  further  support  the  col¬ 
lective  interaction  concept. 

The  mechanism  of  noise  generation  is  also  investigat¬ 
ed  using  far  field  measurements  at  90*to  the  jet  axis. 
The  measured  far  field  power  spectra  exhibit  both  high  and 
low  frequency  peaks.  These  peaks  are  the  signature  of  a 
well  organized  flow  field  between  the  nozzle  and  the 
plate.  The  high  frequency  noise  sources  are  confined  to  a 
short  distance  from  the  nozzle  exit  (due  to  the  effect  of 
the  collective  interaction),  and  radiate  sound  directly  to 
the  far  field.  The  enhanced  large  scale  structures  (due 
to  the  feedback  and  vortex  stretching!  radiate  low  fre¬ 
quency  noise  when  they  impinge  on  tie  plate. 
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CHAPTER  1 


INTRODUCTION 


The  impinging  of  a  turbulent  jet  on  a  surface  is  an 
interesting  fluid  mechanics  research  subject.  The  problem 
has  a  diverse  flow  field  that  includes  a  free  turbulent 
jet  before  impingement,  augmented  vortices  interacting 
with  the  surface,  and  a  developed  wall  jet.  Furthermore, 
at  high  subsonic  impinging  jet  speeds,  the  flow  develops 
self-sustained  oscillations  accompanied  by  an  intense 
screech  tone.  The  main  motivation  of  the  present  study  is 
to  arrive  at  a  better  understanding  of  the  feedback  me¬ 
chanism  responsible  for  the  resonance  phenomenon  at  high 
jet  speeds. 

The  jet  impingement  process  is  also  important  in  many 
applications.  For  example,  short  take-off  and  landing 
(STOL)  aircraft  utilizes  an  externally  blown  flap  techni¬ 
que.  Hot  gasses  that  exhaust  from  the  engine  at  high 
speeds  are  deflected  by  direct  impingement  on  the  flaps  to 


* 
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create  extra  lift  during  take-off.  Problems  encountered 
In  this  design  Include  heat  transfer  and  fatigue  due  to 
excessive  loading  on  the  flaps,  and  high  levels  of  noise 
radiation . 


1.1  REVIEW  OF  PREVIOUS  INVESTIGATIONS 


According  to  Refs.  1  to  4,  the  flow  field  of  an  Imp¬ 
inging  jet  Is  divided  into  three  regions.  A  free  jet  re¬ 
gion,  where  the  mean  jet  characteristics  remain  unchanged 
despite  the  presence  of  the  plate,  an  Impinging  region, 
and  a  developed  wall  jet  region. 


Donalson,  et  al  .  (1971)  scaled  the  heat  transfer  and 
pressure  gradient  at  the  impinging  region  with  the  mean 
properties  of  the  free  jet  at  the  impinging  plane,  as  if 
the  plate  was  removed.  Giralt,  et  al .  (1977)  preferred 
to  use  the  free  Jet  properties  in  a  plane  upstream  of  the 
plate  to  scale  the  impinging  region  properties.  The  plane 
was  chosen  such  that  the  centerline  velocity  is  98%  of  its 
free  value.  In  this  way,  they  accounted  for  the  change  in 
entrainment  at  the  region  of  strong  interaction  near  the 
plate,  as  compared  to  the  scaling  of  Donalson,  et  al . 
(1971).  Preisser  and  Block  (1976)  noted  that  the  surface 
pressure  fluctuations  normalized  with  the  dynamic  head  in 
the  wall  jet  region  are  independent  of  the  la  oh  number  and 
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the  distance  from  the  nozzle  to  the  plate.  Gutmark,  et 
•1.  (1976)  measured  velocity  fluctuations  in  a  plane  jet 
near  its  impingement  on  a  flat  plate.  Spectral  measure¬ 
ments  showed  that  the  plate  plays  the  role  of  a  low-pass 
filter  for  the  turbulent  energy.  A  neutral  scale  was  de¬ 
tected  at  which  turbulent  energy  is  neither  augmented  nor 
attenuated.  At  higher  frequencies  the  turbulence  is  at¬ 
tenuated  due  to  viscous  dissipation,  while  at  lower  fre¬ 
quencies  the  turbulence  is  augmented  due  to  vortex 
stretching  as  the  plate  is  approached.  The  convection 
velocity  of  turbulent  vortical  eddies  on  the  plate  was 
measured  by  Strong,  et  al .  (1  967)  via  two-point 
cross-correlation.  The  velocity  of  convected  eddies  was 
found  to  vary  along  the  plate. 

For  high  subsonic  impinging  jets,  an  audible  screech 
tone  was  reported  in  Refs.  7,  8,  and  9.  Westley,  et  al . 
(1972)  observed  the  screech  tone  and  speculated  a  similar 
mechanism  to  that  of  an  under-expanded  supersonic  jet. 
Wagner  (1971),  and  later  Neuworth  (1973),  reported  that 
the  screech  tone  occurs  at  jet  Mach  numbers  greater  than 
.6  when  the  plate  is  less  than  six  diameters  from  the  noz¬ 
zle  . 


Supported  by  visualizations,  Wagner  (1971)  and  Neu¬ 
worth  (1973)  adopted  the  following  model  to  explain  the 
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feedback  mec’-anism  responsible  for  thie  resonance  phenome¬ 
non.  The  observed  standing  wave  In  the  jet  core  was  as¬ 
sumed  to  result  from  the  superposition  of  a  downstream 
convected  wave  due  to  the  large  turbulent  eddies,  and  an 
upstream  propagating  wave  emanating  from  the  stagnation 
region  of  the  plate.  The  jet  core  simulated  an  acoustic 
wave-guide,  with  the  shear  layer  as  its  boundaries.  The 
upstream  acoustic  wave  was  assumed  to  propagate  in  a  di¬ 
rection  inclined  to  the  jet  axis,  and  to  reflect  from  the 
shear  layer  with  phase  reversal.  The  upstream  waves  re¬ 
ached  the  nozzle  and  forced  the  shear  layer  there  to  form 
a  new  series  of  eddies  that  convected  downstream.  A  mini¬ 
mum  Jet  Mach  number  for  the  resonance  to  occur  was  found, 
at  which  the  shear  layer  acted  as  a  sound  hard  boundary, 
which  means  the  forcing  of  the  shear  layer  near  the  nozzle 
was  not  possible.  Also,  a  minimum  limit  for  the  resonance 
frequency  was  derived  below  which,  the  acoustic  waves  were 
so  inclined  to  the  Jet  axis,  that  no  upstream  propagation 
towards  the  nozzle  took  place.  Irrespective  of  the  Mach 
number  or  the  frequency,  resonance  disappeared  when  the 
plate  was  one  to  two  diameters  beyond  the  potential  core 
of  the  jet.  The  reason  for  this,  was  the  disintegration 
of  ring  vortices  at  these  downstream  locations. 

Wagner  (1971)  and  Neuworth  (1973)  showed  that  the  re¬ 
sonance  frequency  decreases  as  the  plate  is  moved  down- 


stream  until  It  reaches  a  minimum  and  Jumps  to  a  higher 
value.  Vagner  predicted  that  the  decrease  in  frequency  is 
due  to  an  increase  in  the  spacing  between  the  vortices, 
and  the  frequency  jump  occurs  due  to  the  addition  of  a  new 
vortex  in  the  flow  field  between  the  nozzle  and  the  plate. 
In  a  similar  fashion,  Neuworth  speculated  that  the  sum  of 
the  individual  periods  of  the  vortices  and  the  upstream 
waves  is  a  number  that  r  eras  ins  unchanged  between  consecu¬ 
tive  jumps.  The  jump  takes  place  as  a  result  of  an  incre¬ 
ase  by  one  of  the  total  number  of  periods  in  the  loop. 

The  resonance  frequency  Jumps  and  the  observed  en¬ 
hancement  of  large  coherent  eddies  furnish  a  striking  re¬ 
semblance  to  similar  shear  layer-induced  self-sustained 
oscillations  in  the  flow  over  cavities  and  edge  tones.  A 
close  look  to  theoretical  and  experimental  work  in  these 
areas  would  be  helpful. 

Most  of  the  models  for  flow  over  cavities  and  edge 
tones  seem  to  have  the  following  general  features 
[Rockwell  and  Naudascher  (1979)]:  the  shear  layer  inter¬ 
acts  with  the  solid  boundary  (a  corner  or  a  vedg,e)  down¬ 
stream;  acoustic  or  hydrodynamic  disturbances  are  gener¬ 
ated  and  trigger  the  instability  of  the  shear  layer  up¬ 
stream.  Powell  (1961a)  modeled  the  acoustic  source  at  the 
downstream  boundary  as  a  distribution  of  dipoles.  Later 
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(1961b)  he  showed  that  the  radiated  field  of  dipoles  and 
vorticity  are  equivalent.  The  flow  field  was  presented  as 
a  control  system  with  a  complex  transfer  function.  For 
self-sustained  oscillations,  the  magnitude  of  the  transfer 
function  is  unity,  and  the  phase  angle  is  2N7t *  Rossiter's 
(1964)  empirical  model  suggested  that  the  vortices  inter¬ 
act  with  the  downstream  wall  of  the  cavity  and  the  gener¬ 
ated  acoustic  waves  propagate  upstream  inside  <:he  cavity, 
causing  the  shedding  of  new  vortices  near  the  upstream 
edge.  The  sound  source  in  the  mathematical  model  of  Tam 
and  Block  (1978)  was  chosen  proportional  to  the  fluctua¬ 
tions  of  the  shear  layer  at  the  downstream  corner. 
Consequently,  the  acoustic  waves  are  generated  due  to  ^he 
interaction  of  the  inviscid  flow  with  the  downstream 
corner.  The  model  accounted  for  finite  shear  layer  thick¬ 
ness  effects,  and  acoustic  reflections  from  the  bottom  and 
upstream  cavity  walls  in  calculating  the  cavity  pressure 
oscillations.  Sarohia  (1976)  adopted  a  different  model 
for  the  low  speed  flew  over  cavities.  Here  cavity  oscil¬ 
lations  result  from  propagating  disturbances  which  get  am¬ 
plified  along  the  cavity  shear  layer,  rather  than  from 
acoustic  feedback.  Hot  wire  measurements  and  visualiza¬ 
tions  showed  that  the  shear  layer  deflects  in  and  out  of 
the  cavity  near  the  downstream  corner.  This  lateral  mo¬ 
tion  of  the  shear  layer  resulted  in  a  periodic  shedding  of 
vortices  from  the  downstream  corner  at  the  resonance  fre- 


quency.  The  shedding  of  vortices  at  the  downstream  corner 
was  rather  shown  as  "clipping"  and  "nonclipping"  of  vor¬ 
tices  in  the  pictures  of  Rockwell  and  Knlsely  (1978). 

Flow  visual  1 za t ions  [Rockwell  and  Naudascher  (1  979  )) 
showed  that  the  presence  of  a  downstream  solid  boundary 
tends  to  organize  and  enhance  large  scale  structures  in 
the  turbulent  flow.  This  enhancement  takes  place  mostly 
through  a  feedback,  or  sometimes  just  through  stretching 
of  vortices  due  to  the  diverging  streamlines.  A  study  of 
large  coherent  structures  in  turbulent  shear  flows  is, 
therefore,  essential. 

The  visualizations  of  Brown  and  Roshko  (1971)  showed 
for  the  first  time  that  large  scale  two-dimenpional  vor¬ 
tices  make  up  the  fully  turbulent  two-dimensional  mixing 
layer.  Laufer  (1975)  suggested  that  these  quasi-ordered 
structures  are  deterministic,  in  the  sense  that  they  have 
a  characteristic  shape,  Bize  and  convection  motion  that 
can  be  determined  within  a  relatively  small  standard  devi¬ 
ation.  A  somewhat  different  picture  of  turbulence,  there¬ 
fore,  emerges  that  reflects  the  double- struct  tired  nature 
of  turbulent  flows:  time  dependent  quasi-ordered  large 
scale  structures,  and  fine-scale  random  structures. 
Numerous  experiments  have  been  conducted  to  document  the 
characteristics  of  these  deterministic  orderly  structures. 
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Lau,  et  al.  (1972)  reported  results  from  different 
experimental  investigations  that  consistently  yield  a 
value  of  0.6U  for  the  convection  velocity  of  large  scale 
structures#  They  modeled  the  shear  layer  as  an  axial 
array  of  discrete  vortices  spaced  at  about  1 . 2 5 d  (where  d 
is  the  nozzle  diameter).  Measurements  were  made  in  the 
potential  core,  shear  layer  and  entrainment  region  of  the 
jet.  Their  model  helped  to  explain  the  phase  relation¬ 
ships  between  measured  fluctuating  pressure  and  components 
of  fluctuating  velocity.  Petersen  (1978)  measured  the 
convection  velocity  of  large  scale  structures  from 
cross-correlation  of  near  field  pressure  signals.  The 
convection  velocity  was  found  to  be  constant  and  equal  to 
. 6U  over  the  first  five  diameters.  It  then  decreased  con¬ 
tinuously  further  downstream.  The  spatial  stability  theo¬ 
ry  of  Michalke  (1972)  showed  that  the  wavelength  of  the 
initial  instability  is  determined  by  the  initial  shear 
layer  thickness,  which  depends  on  the  exit  Reynolds 
number.  The  most  unstable  wave  has  a  frequen-y  which  cor¬ 
responds  to  f gj /U  ■  .017  (where  is  the  momentum  thick¬ 
ness)  and  a  phase  velocity  equal  to  .  6U  (same  as  the  value 
measured  for  the  convection  velocity  of  the  large  scale 
structures  downstream).  Crow  and  Champag ne ( 1 97 1  )  observed 
visually  the  instability  waves  of  the  thin  laminar  bounda¬ 
ry  layer  leaving  the  nozzle  lip,  and  orderly  large  scale 


puffs  further  downstream.  The  Instability  waves  scaled 
with  the  thickness  of  the  boundary  layer  [Michal ke ( 1  9 71  )  )  . 
In  contrast,  the  large  puffs  scaled  with  the  notale  diame¬ 
ter,  and  were  forming  at  a  constant  average  Strouhal 
number  fd/U«0.3  .  They  observed  that  the  Jet  was  most 
unstable  when  the  forcing  frequency  corresponded  to  a 
Strouhal  number  of  0.3. 

The  growth  of  the  shear  layer  is  attributed  mainly  to 
"pairing"  between  consecutive  large  vortical  structures  as 
was  first  pointed  out  by  Winant  and  Browand  (1974). 
Laufer  (1974)  hinted  that  turbulent  energy  production 
which  i 8  the  energy  transfer  from  large  structures  to 
small  scale  turbulence,  is  phase  related  to  the  pairing 
between  large  eddies.  Indeed,  that  was  observed  later  in 
the  disintegration  of  large  vortical  rings  and  the  produc¬ 
tion  of  small  scale  turbulence  during  the  pairing,  espe¬ 
cially  near  the  end  of  the  potential  core  of  a  round  jet 
[Browand  and  Laufer  (1975)). 

This  review  is  now  completed  by  a  discussion  of  the 
far  field  noise  of  impinging  jets.  The  far  field  noise, 
as  described  by  Curl'6  equation  [see  Ref.  25),  is  the  sum 
of  a  volume  integral  of  Lighthill's  stress  tensor  over  the 
flow  field,  and  a  surface  integral  of  pressure  fluctua- 
tions  over  the  plate.  Powell  (l06n>  applied  Curl's  equa- 
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tlon  to  both  the  real  flow  system  and  an  Image  system* 
The  result  is  that  the  volume  integral  over  the  image  pro¬ 
vides  the  same  contribution  to  the  far  field  noise  as  does 
the  surface  integral  over  the  actual  surface.  Since  the 
image  system  i6  a  pure  reflection  of  the  real  flow, 
Powell's  model  indicates  that  the  role  of  the  plate  is 
passive,  i.e.,  merely  reflecting  the  sound  generated  by 
the  flow. 

Preisser  and  Block  (1976)  reported  levels  of  measured 
far  field  noise  higher  than  those  pertaining  to  the  free 
jet.  The  Increase  in  the  noise  level  is  a  function  of  the 
angle  and  the  nozzle-plate  distance.  Cross-spectra 
between  the  surface  and  the  far  field  showed  that  the 
outer  portion  of  the  impinging  region  (between  d  and  3d 
from  the  stagnation  point)  is  the  major  contributor  of 
noise.  It  is  worth  mentioning  here  that  the  present  study 
will  show  that  large  coherent  structures  impinge  on  the 
plate  surface  in  a  region  between  one  to  two  nozzle  diame¬ 
ters  from  the  stagnation  point.  Preisser  (1979)  showed 
later  that  far  field  spectra  peaks  at  St^  *  .3.  These 
peaks  have  higher  levels  than  those  in  the  spectra  of  a 
free  jet.  He  suggested  that  the  presence  of  the  plate  en¬ 
hanced  the  noise  generated  by  large  coherent  structures. 

The  experimental  studies  of  Preisser  and  Block  (1976) 
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and  Prelsser  (1979)  could  only  localise  the  impinging  re¬ 
gion  as  a  major  apparent  noise-producing  region.  The 
question  of  whether  the  noise  is  originating  in  a  flow  re¬ 
gion  near  the  surface  and  is  being  reflected  by  the  sur¬ 
face,  or  if  it  is  being  generated  by  the  surface  itself, 
could  not  be  answered.  The  exact  role  of  the  plate  in  the 
noise  production  can  only  be  determined  by  evaluating  the 
volume  integral  over  the  flow  in  Curl's  equation.  Such  a 
calculation  requires  near  field  as  well  as  surface  meas¬ 
urements  [Pan  (1975)]. 


1.2  MOTIVATION  AND  SCOPE  OF  THE  PRESENT  STUDY 


The  first  objective  of  the  present  investigation  was 
to  study  the  resonance  phenomenon  that  occurs  in  a  high 
subsonic  impinging  Jet.  The  goal  was  to  explore  a  possi¬ 
ble  feedback  mechanism  between  the  nozzle  and  the  plate, 
and  to  determine  the  role  of  large  coherent  structures  in 
that  mechanism.  The  second  objective  was  to  study  the  far 
field  noise  of  the  impinging  jet  in  general,  i.e.,for  both 
resonant  and  nonresonant  jets  (covering  a  range  of 
.3<M<.9,  and  3<Xg/d<7.5).  Our  goal  was  to  locate  major 
noise-producing  regions,  and  to  study  the  noise-generation 
mechani sms . 

The  present  work  is  based  on  pressure  measurements 
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taken  in  the  far  field,  near  field,  and  on  the  surface  of 
the  plate.  Fig.  (1-1)  ahovs  the  coordinate  system  used 
for  the  jet-plate  configuration  along  with  a  schematic  di¬ 
agram  of  the  pressure  transducers.  Visualizations  in 
water  were  used  to  study  the  general  behaviour  of  coherent 
vortical  structures  during  impingement,  and  are  presented 
in  the  first  section  of  Chapter  3.  The  rest  of  the  inves¬ 
tigation  deals  with  measurements  in  air.  Data  pertaining 
to  pressure  loading  and  eddy  convection  velocity  on  the 
plate  are  included  in  Chapter  3.  Conclusive  experimental 
evidence  for  the  existence  of  convected  large  scale  struc¬ 
tures,  and  upstream  acoustic  waves  generated  from  the  im¬ 
pingement  of  these  structures  on  the  plate  is  presented  in 
Chapter  4.  Also,  a  simple  mathematical  model  purporting 
to  explain  the  particular  behaviour  of  measured  near  field 
correlations  is  included  in  Chapter  4.  A  phase  lock  near 
the  nozzle  exit  between  the  two  waves  necessary  for 
self-sustained  flow  oscillations  is  established  in  Chapter 
5.  A  new  phenomenon  named  "collective  interaction"  is 
also  described  in  connection  with  the  forcing  of  the  shear 
layer  near  the  nozzle  by  the  upstream  acoustic  waves.  The 
chapter  is  concluded  by  a  mathematical  analysis  of  the 
respective  efficiencies  of  acoustic  forcing,  when  the 
shear  layer  is  excited  from  the  ambient  region  and  from 
the  jet  core  region. 
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The  role  of  coherent  flow  including  both  the 
high-frequency  sheer  layer  instability  waves  near  the  not- 
tle  and  the  downstream  impinging  large  scale  structures, 
in  the  noise  generation  is  discussed  in  Chapter  6. 


CHAPTER  2 


EXPERIMENTAL  FACILITIES 

The  experimental  facilities  and  data  processing  meth¬ 
ods  used  in  the  present  investigation  are  discussed  in 
this  chapter.  A  description  of  the  low  speed  impinging 
water  jet  and  associated  visualisation  techniques  are  pre¬ 
sented  first.  Then  the  high  speed  air  jet  facility  where 
most  of  the  measurements  were  taken  will  be  described. 
The  analog  and  digital  data  processing  will  be  discussed 
in  the  last  section. 

2.1  WATER  JET 

The  water  jet  set  up  is  shown  in  Fig.  (2-la)  and 
schematically  in  Fig.  (2-lb).  The  nozzle  assembly  con¬ 
sisted  of  a  stilling  chamber  and  a  detachable  nozzle.  The 
flow  enters  the  stilling  chamber  through  a  conical 
diffuser  with  three  screens,  then  passes  through  a  honey¬ 
comb  straightner  to  the  nozzle.  The  axisymmetric  nozzle 
has  a  contraction  ratio  cf  1.:.'  and  an  exit  diameter  of 


pi 
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1.5  In.  The  jet  exit  velocities  ranged  from  5  in/sec.  to 
15  in/aec.  with  the  corresponding  Reynolds  numbers,  Re  « 

,  of  5000  to  15000.  The  turbulence  level  was  measured 
on  the  jet  centerline  just  upstream  of  the  nozzle  exit  and 
is  reported  in  Ref.  24.  The  rms  longitudinal  fluctuation 
level  is  about  .5X  at  a  jet  speed  of  5  in/sec.,  and  is 
lower  at  higher  jet  speeds. 

In  operation,  the  jet  is  immersed  vertically  in  a  32 
x  32  x  48  in.  plexiglas  tank.  A  centrifugal  pump  capable 
of  providing  0.55  gal/sec.  against  a  head  of  30  feet  of 
water  is  used  to  circulate  the  water  in  this  closed  sys¬ 
tem.  The  jet  impinges  on  a  23  x  23  x  1/2  in.  plexiglas 
plate.  The  ratio  of  the  surface  aren  of  the  plate  to  the 
tank's  cross-sectional  area  is  approximately  1:2.  The  un¬ 
iform  diffusion  of  the  dye  over  the  whole  tank,  above  and 
under  the  plate,  indicated  no  significant  plate  blockage 
effect  on  the  jet  flow.  The  plate  is  supported  from  its 
corners  by  four  1.0  in.  PVC  rods,  which  are  connected  to 
a  wooden  frame  on  top  of  the  tank.  The  wooden  frame  is 
also  used  to  support  the  nozzle  assembly.  The  nozzle  to 
plate  distance  is  changed  by  lowering  the  plate  through 
1/2  in.  steps.  After  impingement,  the  flow  is  collected 
in  an  exhaust  diffuser  at  the  bottom  of  the  tank.  The 
flow  rate  is  monitored  by  a  Fisher-Porter  flow  meter  with 
a  gate  valve  providing  the  required  throttling. 
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The  flow  was  visualised  in  two  ways: 

1 .  Colored  shear  layer 

Dye  was  introduced  into  the  nozzle  wall  boundary 
layer  through  a  series  of  fine  slots.  The  dye  used  was 
common  red  food  coloring.  This  is  a  convenient  coloring 
agent  since  it  can  be  cleared  by  ordinary  laundry  bleach. 
Still  photographs  and  films  were  taken.  Back-lighting 
through  a  transluscent  screen  was  uBed,  and  was  found  to 
give  the  best  lighting  conditions. 

2 •  Hydrogen  bubble  technique 

Hydrogen  bubbles  were  generated  from  a  thin  platinum 
wire  that  was  stretched  across  the  exit  plane  of  the  noz¬ 
zle.  The  wire  was  a  cathode  in  an  electric  circuit,  and 
was  connected  to  a  pulse  generator  and  a  power  amplifier. 
The  anode  was  a  long  rod  immersed  in  the  water  near  the 
wire.  An  electrolyte,  table  salt,  was  added  to  the  water 
in  a  quantity  selected  to  maximize  the  production  of  hy¬ 
drogen  bubbles.  A  .002  in.  diameter  platinum  wire  and 
voltage  pulses  of  approximately  40  volts  in  amplitude  were 
used  to  produce  a  sequence  of  hydrogen  bubble  lines.  The 
wire  was  small  enough  to  produce  the  desired  bubble  size 
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and  yet  strong  enough  to  withstand  drag  forces  due  to  the 
jet  flow.  The  axial  slice  of  the  jet  was  illuminated 

1 

j 

using  either  a  commercial  spot  lamp  with  a  vertical  slit^ 
or  a  stroboscope  lamp  and  a  system  of  lenses  with  a  verti- 1 
cal  slit.  The  stroboscope  was  triggered  by  the  pulse  gen¬ 
erator,  and  this  setup  was  used  for  multi-exposure  still 
pictures.  The  room  was  completely  darkened  to  provide 
good  contrast  between  the  bubbles  and  the  background.  The 
level  of  illumination  was  minimal  for  photography.  The 
negatives  of  both  the  still  pictures  and  movie  films  had 
to  be  pushed  to  as  high  as  ASA  3000  in  the  developing  pro¬ 
cess. 


2.2  AIR  JET 


2.2.1  Anechoic  Chamber 

The  anechoic  chamber  at  the  University  of  Southern 
California  was  constructed  by  Eckel  Industries.  The  inte¬ 
rior  walls  are  faced  with  fiberglas  wedges  which  have  a 
depth  of  16  1/2  in.,  resulting  in  a  lower  cut-off  frequen¬ 
cy  of  150  Hz.  The  interior  measures  15'  8"  vide  by  20*  4" 
long  and  by  11*  11"  high.  Eight  inches  above  the  floor 
wedges  a  tensioned  cable  is  suspended  from  the  walls  to 
permit  easy  access  to  the  chamber.  The  jet  flow  is  vented 
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to  the  outside  by  a  4*  x  4*  opening  with  a  baffled  cut.  At 
various  convenient  locations  in  the  chamber  threaded  cou-  | 
plings  for  standard  1  1/2  in.  I.P.S.  pipes  are  mounted 
to  provide  a  base  for  mounting  the  flat  plate  and  its  tra- 

£ 

verse  mechanism.  The  chamber  is  mounted  on  large  coil 
springs  to  assure  a  complete  isolation  from  any  vibra¬ 
tions.  A  schematic  of  the  chamber  and  jet  facilities  is 
shown  in  Fig.  (2-2). 

2.2.2  Jet  Facilities 

i 

The  air  supply  system  consists  of  a  2500  psi,  142 
3 

j-  ft.  /min.,  Chicago  pneumatic  air  compressor  which  sup¬ 

plies  air  to  five  storage  pressure  vessels.  The  compres¬ 
sor  has  two  mechanical  filters  and  a  chemical  dehydration 
,  system  which  provides  dry  air  with  a  dew  point  of  -4  8  F. 

The  compressor  is  used  to  fill  the  storage  tanks  and  is 
not  operated  during  test  runs. 

The  air  flow  from  the  pressure  vessels  is  controlled 
by  a  low-noise,  self-drag  valve  built  by  Control  Compo¬ 
nents  Inc.  of  Irvine,  California.  The  valve  is  opened  by 
a  Moore  3-15  psi  positioner  which  receives  a  pneumatic 
pressure  signal  from  a  Honeywell  Batch-O-Line  Controller. 
The  controller  was  operated  manually,  while  the  settling 

chamber  pressure  time  history  was  recorded  on  a  strip 
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chart.  A  Gould  Brush  Model  220  dual  channel  chart  record¬ 
er  was  used,  which  was  connected  to  a  Valldyne  Model  DP 
15TL  pressure  transducer.  The  variation  of  the  stagnation 
pressure  during  the  measurements  was  maintained  to  within 
22  of  the  mean  guage  pressure.  The  above  setup  provided 
excellent  monitoring  of  the  mean  flow  pressure  and  this 
proved  to  be  critical  for  the  present  investigation:  we 
observed,  even  audibly,  that  relatively  small  fluctuations 
in  the  stagnation  pressure  in  the  order  of  6  to  92  can 
significantly  alter  the  resonance  at  high  subsonic  jet 
speeds . 

After  discharging  from  the  control  valve,  the  air 
flows  through  a  500  KW  Hynes  Radi-Fin  electric  heater. 
The  heater  was  not  used  in  the  present  investigation. 
Immediately  after  the  heater  the  flow  passes  to  the  set¬ 
tling  chamber  through  a  diffuser  that  is  fitted  with 
screens  to  prevent  flow  separation.  To  reduce  turbulence 
five  screens  and  a  4  in.  layer  of  steel  wool  are  in¬ 
stalled  inside  the  16  15/16"  diameter  settling  chamber. 
The  turbulence  level  downstream  of  the  last  screen  has 
been  measured  to  be  0.42.  This  is,  of  course,  further  re¬ 
duced  by  the  contraction  downstream  of  the  settling 
chamber • 

For  the  present  study  a  subsonic  axisymmetric  nozzle 


of  1"  diameter  was  used,  yielding  a  contraction  ratio  of 
approximately  17:1.  During  this  investigation  the  jet 
Mach  number  was  restricted  to  the  range  0.3  to  0.9.  More  * 
details  about  the  Jet  facilities,  the  contraction  section  1 
downstream  of  the  settling  chamber,  and  the  nozzle  are  re¬ 
ported  in  Ref.  29. 


2.2.3  The  Flat  Plate  and  the  Supporting  Mechanism 

A  flat  steel  plate  (30"x30"x . 25 ")  provided  a  large 
surface  for  the  impinging  jet.  Fig. (2-3)  shows  a  picture 
of  the  nozzle  and  the  plate  inside  the  anechoic  chamber. 
First  holes  were  drilled  in  the  plate  at  dynamic  and  stat¬ 
ic  pressure  measurement  sites.  The  details  of  these  holes 
and  their  approximate  locations  are  shown  in  Fig .( 2-4 ) .The 
exact  locations  of  these  holes  are  shown  in  Table-1.  The 
plates  surface  facing  the  jet  was  then  ground  to  «  commer¬ 
cial  mirror  finish  (16  raicroinch  surface  finish).  Finally 
the  surface  was  treated  with  black  chrome  for  corrosion 
resistance.  The  unused  dynamic  pressure  holes  in  the 
plate  were  filled  by  dummy  transducers. 

The  holes  of  the  pressure  taps  had  clean  sharp  edges 
at  the  surface.  With  both  pressure  and  dummy  transducers 
flush  mounted  on  the  plate,  it  was  believed  that  the  re¬ 
sulting  smooth  surface  would  r.o  t  interfer  with  the  flow 
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near  the  plate 


The  plate  was  mounted  on  a  traverse  mechanism  as 
shown  In  Fig. (2-5).  The  vertical  motion  was  provided  by  a 
motor  that  slides  the  plate  over  a  Thomson  Linear  Ball  Be¬ 
aring  (Type  SPB-12-OPN)  with  shaft  assemblies  Type 

SA12-30.  The  plate  could  rotate  in  a  pitching  mode, and  be 

•  • 

fixed  at  any  angle  to  the  jet  axis  between  90  to  180.  In 
the  present  study,  normal  jet  impingement  only  was  used. 
The  plate  was  connected  by  a  simple  truss  to  a  base.  The 
base  was  connected  to  eight  1.5  in.  I.P.S.  pipes  that 
were  screwed  tightly  into  the  floor  of  the  chamber.  The 
base  allowed  manual  horizontal  transverse  motion  and  mo¬ 
torized  axial  motion.  The  axl8l  motion  was  used  to  vary 
the  nozzle-to-plate  distance  x^.  The  plate  can  move  from 
Xq^O  to  Xq*25  in.  The  vertical  and  axial  motions  were  re¬ 
motely  controlled. 

The  plate  and  its  supporting  mechanism  were  examined 
for  structural  vibrations  during  jet  operation.  The  am¬ 
plitude  of  the  vibrations  of  the  plate  was  measured  by  an 
accelerometer(Bruel&Kjaer  4324).  With  the  highest  speed 
jet,  M*.9  ,the  amplitude  of  vibrations  was  only  1.12x10 
in . 
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2.3  DATA  PROCESSING 


2.3.1  Measuring  and  Recording  Equipment 
Dynamic  pressure  pickups 

The  surface  pressure  fluctuations  were  measured  by 
1/8  in.  Kulite  XTE-1-190-25  pie2oresistive  pressure 
transducers  with  type  B  screens.  From  previous  experience 
difficulties  were  anticipated  during  installation  and  re¬ 
moval  of  the  transducers.  The  screen  was  glued  to  the 
transducer's  tip  and  frequently  would  be  pried  during  re¬ 
moval  because  it  has  the  same  diameter  as  the  transducer 
tip  and  precise  alignment  was  difficult.  A  pressure  port 
for  the  transducer  was  designed.  The  transducer  and  its 
port  [Fig .  (2-6)]  were  always  used  as  a  single  piece  dur¬ 
ing  the  mounting  and  removal  from  the  plate.  In  this  way 
the  transducers  were  protected  from  possible  damage. 


The  natural  frequency  of  the  Kulite  transducer  is  in 
excess  of  100  KHz.  A  preamplifier  circuit  was  built  for 
each  Kulite  and  is  shown  in  Fig.  (2-7).  The  first  stage 
of  the  circuit  is  a  voltage  follower  to  reduce  the  output 
impedance  of  the  transducer.  The  second  stage  is  an  am¬ 
plifier  which  provides  a  gain  of  200,  500  or  1000. 
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Near  field  microphones 

Four  1/8"  B  &  K  Type  4138  condenser  microphones  were 
used.  They  had  a  frequency  response  from  7  Hz  to  1 A 0  KHz 
and  a  range  of  76-185  dB. 

Far  field  microphones 

A  1/4"  B  &  K  Type  4135  condenser  microphone  was  used. 
It  had  a  frequency  response  from  4  Hz  to  100  KHz  and  a 
range  of  59-1 64  dB . 

Microphone  preamplif iers 

Five  B  &  K  Type  2618  preamplifiers  were  used  with  the 
near  field  and  far  field  microphones.  The  frequency  res 
ponse  of  these  preamplifiers  were  2  Hz  -  200  KHz  and  could 
provide  a  gain  of  -20,  0,  or  +20  dB. 

Microphone  power  supply 

A  7-  channel  B  &  K  Model  222  power  supply  was  used. 
A  gain  of  up  to  +40  dB  can  be  provided  by  the  power  supply 
for  each  channel . 
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Microphone  calibrator 

B  &  K  Type  4220  plstonphone  calibrator  was  used  to 
calibrate  the  B  &  K  microphones. 

Magnetic  tape  recorder 

During  each  run,  the  data  were  recorded  on  a  14- 
channel  Hewlett-Packard  Model  3955A  type  recorder. 
Measurements  of  surface  pressure  fluctuations  were  record¬ 
ed  in  FM  mode  with  frequency  bandwidth  from  DC  to  20  KHz 
at  a  tape  speed  of  60  ips.  The  near-field  and  far-field 
data  were  recorded  in  direct  mode  with  response  from  300 
Hz  to  300  KHz  at  the  same  tape  speed. 

On-line  equipment 

During  and  after  the  run  the  signals  were  displayed 
on  the  screen  of  a  two  channel  Hewlett  Packard  Type  121 7A 
oscilloscope.  The  rras  of  the  pressure  signals  were  meas¬ 
ured  from  a  Hewlett-Packard  400E  AC  Voltmeter.  A 
Hewlett-Packard  3721A  correlator  was  used  to  give  on-line 
autocorrelation  and  cross-correlation  of  the  recorded  sig¬ 
nals.  The  delay  offset  provided  in  the  correlator  was 
helpful  for  observing  long  optimum  time  delays  in  the 
cross-correlation  between  near-field  and  far-field  sig- 
ials . 
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2.3.2  Digital  Data  Procerjing 
Digitization 

Most  of  the  analysis  in  the  present  study  was  done 
digitally.  The  analog  data  were  digitized  and  processed 
on  a  Digital  Equipment  Corporation  PDP-11/55  minicomputer. 
The  subsystem  used  for  digitization  consists  of  16  sample 
and  hold  circuits  and  an  analog  to  digital  converter  in¬ 
terfaced  with  the  minicomputer.  The  maximum  sampling  rate 
attained  in  this  study  was 

h  ■  13  x  (number  of  channels)  +  20  microseconds 

The  maximum  sampling  rate  was  found  to  depend  on  many  par¬ 
ameters  including  how  busy  the  computer  was  during  digiti¬ 
zation,  the  block  size  of  digitized  data  written  to  mag¬ 
netic  tape,  and  the  quality  of  the  digital  tape.  However, 
an  effective  faster  sampling  rate  could  be  obtained  by 
slowing  down  the  analog  tape  during  the  digitization  pro¬ 
cess.  A  slow  factor  (SLF)  as  much  as  32  was  used  for  some 
cases,  by  using  an  analog  playback  tape  speed  of  1  7/8  ips 
during  digitization.  As  mentioned  before  recordings  were 
performed  at  60  ips  s per c . 
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The  high  frequency  cutoff  from  digitization  [Ref. 
30]  is 

.  (sir)  J_  . 

where  f c  is  the  Nyquist  folding  frequency.  The  frequency 
resolution  is 

B  «  (slf) — L_ —  ,  (a.o 

Nr  .  h 

where  NT  is  the  number  of  points  per  channel  in  the  record 
being  processed.  The  equations  above  show  a  trade  off 
between  the  frequency  range  and  the  frequency  resolution. 
(SLF)  ■  32  was  used  to  look  at  the  high  frequency  compo¬ 

nents  of  a  signal,  such  as  in  measuring  the  shear  layer 
initial  instability  frequency  in  Ch .  5  and  the 

high-frequency  far  field  noise  in  Ch .  6.  Otherwise,  slow 

factors  of  1  and  8  were  used. 

Data  analysis  procedures 

Figure  (2-8)  shows  a  block  diagram  depicting  most  of 
the  digital  data  analysis  procedures  used  in  the  present 
study.  The  limitations  on  the  available  computer  storage 
required  that  the  Fast  Fourier  Transform  (FFT)  and  the 
evaluation  of  the  power  spectrum  components  be  carried  out 
over  one  block  of  digital  data  and  then  be  summed  and 
averaged  over  the  total  number  of  blocks  digitized.  This 
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is,  In  effect  a  high  frequency  filter  with  a  low  frequency 
cutoff  equal  to  the  frequency  resolution  B  given  by  Eqn . 
(2-1). 


The  spectra  and  correlation  functions 

Most  of  the  following  analysis  is  based  on  defini¬ 
tions  and  notations  given  in  Ref.  30.  The  cross  spectrum 
function  for  two  time-signals  i(t)  and  j(t)  is  by  defini¬ 
tion 

%  *  c,j  lf>  -  i  >  <•*-•*«> 

where 

T 

C:  (f)s  '  Be)  dt  ,  U-Jb) 

TBe  ° 

which  is  called  the  coincident-spectral  density  function 
and 

T 

o*  JL- J  &«)  dt  ,  (2.20 

J  Toe  ° 

which  is  called  the  quadrature-spectral  density  function. 
In  Eqn.  (2-2)  i(t,f#B  )  and  j  (t,f,B  )  are  the  signals 
passed  through  a  narrow  band-pass  filter  with  bandwidth  Bg 
and  center  frequency  f;  j* ( t  ,  f , B  )  is  90-degree  phase 
shifted  from  j(t,f,Bg).  T  is  the  record  time  of  both  i(t) 
and  j(t)  and  is  equal  to  1/B. 
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A  more  useful  form  of  Eqn •  (2-2)  is 


1  2  ft 

[c;.«)+  Q..tnl 


c%)  I J  it) 


'•J  l  C:..  »i  J  J 


(2-5) 


where  G.  .(f)  is  a  normalised  spectrum  satisfying 

1  I  J 

J  6.  (f)  df  s  I 

*  <»j 

and  *'?.  .(f)  is  the  phase  angle  function. 

1  9  J 

The  cross-correlation  function  is  given  in  normalized 


form  by 


R  CT)s  Ut)  jtt  +  T) 
'•J  "  r= —  r~ 


(2~4 ) 


I A  +t; 


where 


at)  j(t+T)  =  J  ut)  ju+T)  at 

T  o 

Eqn.  (2-4)  is  evaluated,  as  shown  in  Fig.  (2-8),  by  an 
Inverse  Fourier  Transform  of  the  cross-spectrum  G.  .(f), 

*  9  J 

i  .e . 


R.  .IT) 


nr*  T+Qlf)  s im  JtnfTj  d  T  . 

•i 


If  the  input  signals  in  Fig.  (2-8)  are  the  same,  say 

both  are  i(t),  Eqn .  (2-2)  gives  the  power  spectrum  G^f), 

and  Eqn.  (2-4)  the  autocorrelation  function  R,  . (T) 

1  *  1 


G;(f)  =  -JL  -1-  1  Bt>  dt  (J.S) 

l*t)  TBC  0 


Ri,i  <T>  • 


T 

-i-  j  at)  at*T)  at  . 
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A  rough  estimate  of  the  spectrum  can  be  evaluated 
from  Eqns.  (2-2)  and  (2-3)  for  Bg  *  B  and  is  referred  to 
in  the  text  as  G(f).  If  the  signals  were  of  a  limited 
band-width  white  noise,  the  standard  normalized  error  of 
the  spectral  estimates  would  be  *  1.0,  which  means  that 
the  error  is  of  the  same  order  as  the  estimate  itself. 
However,  most  of  the  signals  measured  in  the  present  study 
are  of  a  narrow-band  random  noise  type  and  is  expected 
to  be  less  than  unity.  This  proved  to  be  true  when  G(f) 
was  compared  with  the  smooth  spectrum  estimates  in  the 
form  of  the  coherence  function  .(f),  which  is  defined 

1  f  J 
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as 

2 

\  .  to  =  __  _  . 

,J  G;  (0  Gj  « ) 

The  barred  quantities  are  the  smooth  spectral  estimates 

averaged  over  frequency  band  of  width  Bg  ■  nB.  The 

changes  in  spectral  estimates  for  different  values  of  n 

were  studied  for  different  samples  of  signals  and  were 

found  to  be  negligibly  small  for  n^.8  .  The  value  for  n  *  8 

was,  therefore,  chosen  in  calculating  .(f)  in  the  pre- 

1 » J 

sent  study. 

The  prewhitening  technique 

As  we  shall  see  later,  at  high  subsonic  jet  speeds 
the  pressure  signals  take  the  form  of  distorted  sine 
waves.  The  spectra  then  are  dominated  by  a  peak  at  the 
resonance  frequency,  while  the  cross-correlation  function 
would  show  a  sinusoidal  behavior  with  no  distinguishable 
optimum  peak.  The  cross-correlation  functions  in  these 
cases  do  not  provide  much  information.  A  data  processing 
method  called  prewhitening  is  quite  useful  for  analyzing 
small  random  signals  superimposed  on  a  high  amplitude  pure 
tone.  A  conventional  prewhitening  technique  [Williams  and 
Purdy  (1970))  separates  the  random  signals  from  the  pure 
tone  by  subtracting  an  in-phase  sine  wave  from  the  origi- 
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nal  signal  by  analog  methods*  A  prewhitening  technique 
vaa  developed  here  using  digital  methods,  and  was  found  to 
be  much  simpler.  A  c ro ss- cor rel a tlon  function  of  two 
prewhitened  signals  was  obtained  by  the  following  sequence 
of  operations: 


1.  Compute  the  individual  power  spectra  G^(f)  and 

G .(f)  and  also  the  cross-spectrum  G.  .(f). 

J  1  »  J 

2.  Remove  the  pure  tone  from  the  above  three  spectra 
by  replacing  the  high  spectral  peak  with  the  averaged 
value  of  the  two  neighboring  frequency  components.  This 
results  in  a  prewhitened  G^(f),  Gj(f)  and  Gi  j(f)» 

3.  Inverse  Fourier  Transform  G^(f)  and  G^(f)  to  get 
the  autocorrelations  R.  ,  (T)  and  R,  .  (r).  Note  that 

^  ^  J  I  J 

R.  At)  and  R,  . (T)  are  the  mean  square  values  of  the  ran- 
I  *  I  J  *  J 

dom  signals  remaining  from  i(t)  and  j(t). 


4.  Inverse  Fourier  Transform  G.  .(f)  and  normalize 

1  9  J 

to  obtain 


R.  (r)  = 

l»j 


i riU(o> 


(2-7) 


Results  derived  from  ~>r  evhi  tening  method  are  shown  in 
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CHAPTER  3 

GENERAL  CHARACTERISTICS  OF  THE  IMPINGING  JET 


A  general  idea  about  the  common  properties  of  an  imp¬ 
inging  Jet  is  presented  here  before  discussing  specific 
details  of  the  feedback  mechanism  and  the  noise  radiation. 
In  the  first  part,  visualizations  in  water  are  used  to 
describe  the  general  behavior  of  large  vortical  rings  dur¬ 
ing  their  impingement  on  a  flat  plate.  The  water  jet  is 
limited  to  a  low  Reynolds  number  range,  Re  ■  5000  to 
15000.  In  the  second  part,  surface  pressure  loading,  eddy 
convection  velocity  on  the  plate  ,and  the  stability  modes 
of  the  impinging  jet  are  discussed  for  air  Je^s  with  a  Re¬ 
ynolds  number  range  of  1.8  x  10^  to  5.3  x  10^. 

3.1  VISUALIZATION  OF  JET  IMPINGEMENT -WATER  JET 

As  mentioned  in  Chapter  J,  the  flow  field  of  the  imp¬ 
inging  jet  is  divided  into  three  regions:  the  free  jet 
region,  the  impinging  region,  and  the  wall  jet  region. 
From  the  dye  and  hydrogen  bubble  techniques,  coherent 


structures  in  the  free  Jet  region  end  the  impinging  region 
were  clearly  visible.  Both  the  dye  and  the  hydrogen  bub¬ 
ble  lines  lost  their  coherent  features  due  to  the  develop¬ 
ing  turbulence  in  the  wall  jet  region. 

The  flow  field  between  the  nozzle  and  the  plate  is 
shown  in  the  dye  picture  of  Fig.  (3-1).  The  flow  is  ini¬ 
tially  laminar  for  approximately  one  diameter  downstream 
of  the  nozzle  exit.  The  shear  layer  rolls  into  axisymme- 
tric  vortical  rings.  These  rings  are  convected  downstream 
until  they  impinge  on  the  plate.  For  the  same  flow  param¬ 
eters,  a  similar  pattern  is  shown  in  Fig.  (3-2)  except 
that  both  the  initial  laminar  column  near  the  nozzle  exit 
and  the  spacing  between  the  vortical  rings  are  larger. 
The  larger  spacing  between  rings  takes  place  through 
’'pairing”  between  two  consecutive  rings.  Observation  of 
the  movies  showed  the  flow  intermittently  jumping  between 
the  two  patterns  of  Fig.  (3-1)  and  (3-2). 

The  two  patterns  described  above  are  further  shown  by 
the  cohernt  flow  motion  presented  by  the  x-t  diagram  of 
Fig. (3-3).  The  diagram  is  constructed  from  a  frame  by 
frame  analysis  of  a  movie.  Each  curve  in  the  figure 
starts  when  a  "bulge”  in  the  initial  laminar  column  is  de¬ 
tected.  These  ''bulges"  evolve  into  vortical  rings  down¬ 
stream.  Both  the  bulge  and  its  corresponding  toroidal 


vortex  is  refered  to  here  as  a  vortical  ring.  Hence, each 
curve  in  the  figure  represents  an  individual  ring  path  or¬ 
iginating  just  downstream  of  the  nozzle  exit  and  terminat¬ 
ing  near  the  plate  before  impinging.  The  of  the  time  axis 
is  arbitrary  and  is  not  indicative  of  the  first  ring  ob¬ 
served  in  the  run.  For  t>2  sec.  pairing  does  not  occur 
in  this  x-t  diagram.  A  vertical  line  at  t"2.4  sec., for 
example,  intersects  ring  paths  at  points  corresponding  to 
the  Instantaneous  locations  of  the  rings.  A  snap  shot  at 
t»2.4  sec.  would  give  a  flow  pattern  similar  to  that  of 
Fig. (3-1).  During  the  period  .5  sec.<t<2  Bee.  pairing 
takes  place.  Here  intersections  between  ring  paths  and 
the  vertical  line  at  t" 1.6  sec.  give  instantaneous  ring 
locations  which  correspond  to  the  pattern  displayed  in 
Fig.  (3-2).  Based  on  the  intersections  between  vertical 
lines  and  ring  paths,  the  following  remarks  can  be  made: 

1.  The  rings  are  approximately  evenly  spaced  between 
the  nozzle  and  the  plate. 

2.  Through  pairing,  the  flow  reorganizes  itself  into 
a  different  pattern,  where  the  spacing  between  consecutive 
rings  is  almost  doubled. 

Browand  and  Lauf er ( 1 975 )  and  later  Pe tersen (  1  978 ) 
have  shown  that  vor.ical  rings  were  qua siper iod ica lly 


spaced  in  a  free  jet.  A  comparison  between  their  x-t  di¬ 
agrams  and  Fig. (3-3)  showed  that  the  rings  seem  to  be  more 
regularly  spaced  for  the  impinging  jet.  Furthermore,  vor¬ 
tices  are  shed  at  a  rathar  constant  frequency  as  one  would 
deduce  from  intersections  between  ring  paths  and  a  hori¬ 
zontal  line  at  x/d<1.0  in  Fig. (3-3).  This  indicates  that 
the  presence  of  the  plate,  especially  near  the  nozzle, 
tends  to  organize  the  flow  field  between  the  nozzle  and 
the  plate.  A  possible  "hydrodynamic"  feedback  could  exist 
at  such  low  subsonic  speeds  [Ho  and  Nosseir  (1  977  ) )  . 
However,  further  investigations  are  needed  before  one 
could  make  a  definite  conclusion  in  that  matter. 

Winant  and  Browand (  1 974  )  were  the  first  to  observe 
pairing  between  vortices  in  the  two-dimensional  shear 
layer.  They  pointed  out  that  the  shear  layer  growth  is 
due  to  pairing  between  these  vortical  structures. 
Petersen(1978)  argued  that  rings  coalesce  when  the  local 
mixing  layer  attains  a  critical  thickness  which  scales 
with  the  wave  length  of  the  vortex  pair.  Pairing  between 
consecutive  rings  is  shown  in  Fig.(3-4)and  is  described  as 
follows:  the  downstream  ring  slows  down  and  grows  slight¬ 
ly  in  diameter,  while  the  upstream  ring  increases  its  vel¬ 
ocity  and  shrinks.  The  upstream  ring  is  then  "swallowed" 
by  the  downstream  one  and  they  both  form  a  larger  ring. 
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The  pairing  process  is  also  shown  in  the  hydrogen 
bubble  picture  of  Fig. (3-5).  The  bubble  lines  near  the 
nozzle  exit  are  equally  spaced  indicating  a  constant  flow 
speed,  and  form  horizontally  straight  lines  indicating  a 
uniform  speed;  these  effects  are  characteristic  of  the 
potential  core.  Each  line  curves  at  the  outer  edge  of  the 
jet  due  to  the  retarded  velocity  in  the  shear  layer.  In 
this  region  the  hydrogen  bubbles  approximately  provide  in¬ 
stantaneous  u-velocity  profiles.  Two  lines  are  shown 
curved  in  opposite  directions  near  the  plate.  The  upper 
curvature  indicates  a  high  convection  velocity  in  the 
"curved"  shenr  layer.  In  fact,  this  is  the  high  velocity 
upstream  ring  as  it  catches  up  with  the  slow  moving  down¬ 
stream  ring  as  mentioned  before.  Browand  and  Laufer(1975) 
described  an  identical  pairing  process  based  on  the  vor¬ 
tices  mean  velocity  from  their  x-t  diagram.  Fig. (3-5) 
further  provides  the  6hape  of  the  velocity  profiles  of  the 
vortices  during  pairing.  In  the  present  experiment  the 
convected  rings  stretch  as  they  impinge  on  the  plate. The 
ring  cores  impinge  at  locations  1.0<  r^/d  <1.5  depending 
on  the  plate  distance  from  the  nozzle.  Fig.  (3-6)  shows 
the  survival  of  vortex  rings  as  they  stretch  to  more  than 
two  nozzle  diameters  after  impingement.  Ho  and  Kovasznay 
(1974)  described  a  similar  process  for  the  impingement  of 
isolated  vortical  rings  on  a  plate.  They  observed  that 
the  ring  radius  increases  rapidly  on  the  plate  to  a  limit- 
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ing  value.  The  vortex  ring  then  disappears  by  dissipa¬ 
tion.  The  limiting  radius  was  found  to  be  a  function  of 
the  downstream  location  of  the  plate.  The  process  is  sim¬ 
ilar  in  the  case  of  the  impinging  jet.  The  difference  is 
that  the  jet  rings  are  superimposed  on  a  mean  flow. 
Therefore,  the  vortical  rings  continue  to  stretch  on  the 
plate  because  of  mean  flow  velocity,  until  they  disinteg¬ 
rate  by  dissipation. 

The  slope  of  ring  paths  of  Fig.  (3-3)  gives  the  con¬ 
vection  velocity  of  these  structures.  The  average  convec¬ 
tion  velocity  of  40  rings  was  found  to  be  .67  U,  which  is 
approximately  the  same  value  measured  by  Browand  and 
Laufer  (1975)  in  the  free  jet.  Furthermore,  it  was  ob¬ 
served  from  the  movie  that  the  velocity  of  ring  cores 
along  the  curved  shear  layer  near  impingement  does  not 
change  significantly.  At  the  mean  time,  the  potential 
flow  in  the  jet  core  decelerates  as  it  approaches  the 
stagnation  point  on  the  plate.  This  is  shown  in  Fig. 
(3-7).  The  distance  between  bubble  lines  decreases  as  the 
flow  approaches  the  plate,  which  indicates  a  decelerating 
flow  near  the  jet  axis.  The  convex  curvature  towards  the 
plate  side  of  the  two  lines  nearest  to  the  plate  indicates 
a  slower  speed  near  the  center  of  the  jet  as  compared  to 
the  outer  shear  layer. 
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The  shear  layer  growth  is  shown  in  the  multiple  expo¬ 
sure  hydrogen  bubble  picture  of  Fig.  (3-8).  The  rate  of  ] 
increase  of  the  shear  layer  thickness  is  faster  near  the  J 
plate.  Note  that  bubbles  produced  by  the  platinum  wire  ■ 
outside  the  Jet  column  in  Fig.  (3-7),  display  nicely  the 
particle  paths  of  the  ambient  entrainment  in  the  multiple 
exposure  of  Fig.  (3-8).  These  entrainment  lines  appear 
clearer  in  the  left  hand  side  of  the  jet  because  of  its 
closer  proximity  to  the  light  source.  '  The  lines  indicate 
that  most  of  the  entrainment  takes  place  near  the  plate, 
especially  at  the  outer  edge  of  the  impinging  region  ,i.e. 
at  the  end  of  the  curved  shear  layer  1.5  <r^/d<  2.0. 
Fairing  between  vortices  on  the  plate  was  observed  to  take 
place  in  this  region  .  The  high  turbulent  energy  produc¬ 
tion  associated  with  the  pairing  [Laufer  (1975)]  is  the 
cause  for  this  high  level  of  entrainment. 

Widnall  and  Sullivan  (1973)  and  Maxworthy  (1976)  ob¬ 
served  azimuthal  instability  of  isolated  vortex  rings. 
The  amplitude  of  the  instability  waves  grew  in  a  direction 
45  to  the  direction  of  ring  propagation.  The  number  of 
waves  was  determined  from  the  condition  of  zero 
self-induced  rotation  of  the  waves  around  the  filament 
axis  in  the  theoretical  model  of  Widnall,  et  al  .  (1  974  ). 

Maxworthy  (1976)  showed  that  the  number  of  instability 
vaves,  which  is  inverse!;,  proportional  to  the  ring  core 
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diameter,  increases  mono tonically  with  the  ring  Reynolds 
number,  Brovand  and  Laufer  (1  9  75)  noticed  enhanced  azimu-  1 
thal  instability  during  the  pairing  of  two  rings  in  the 

i 

free  Jet.  They  argued  that  this  secondary  instability  ex-  1 
plains  the  breakdown  of  the  rings'  coherence  downstream  of 
the  potential  core.  Azimuthal  instability  was  also  ob¬ 
served  in  the  impinging  vortex  ring  and  is  shown  in  Fig. 
(3-9).  However,  the  presence  of  the  plate  and  the  associ¬ 
ated  stretching  of  the  vortices  during  impingement  seems 
to  prevent  or  at  least  delay,  the  breakdown  of  the  coher¬ 
ent  rings  due  to  these  waves. 

At  Re  ^  8000  a  higher  instability  mode  was  observed. 
At  Re  ■  8000,  Fig.  (3-10)  shows  a  spiral  mode  that  was 

observed  to  occur  intermittently  with  the  axisynmetric 
mode.  This  was  also  observed  in  the  free  Jet  by  Browand 
and  Laufer  (1975). 

3.2  PRESSURE  MEASUREMENTS  -  AIR  JET 

The  study  of  the  water  jet  presented  above  was  predo¬ 
minately  qualitative  in  nature.  The  remaining  part  of  the 
present  investigation  deals  with  quantitative  results  der¬ 
ived  from  the  study  of  an  impinging  air  jet.  Surface 
pressure  fluctuations  due  to  the  impingement  of  a  turbu¬ 
lent  jet  are  generally  random  with  respect,  to  time  for  low 


Jet  Mach  numbers,  M.  As  an  example,  the  raw  signal  from  a 
surface  pressure  transducer  at  M  ■  .5  is  shown  in  Fig. 
(3-11).  However,  the  pattern  of  the  signal  changes  com- 

i 

pletely  for  M>.6  and  XQ/d<6.  At  M-0.9  the  signal  can  be 
described  as  a  slightly  distorted  sine  wave  as  shown  in 
the  figure.  This  behaviour  was  usually  accompanied  by  a 
screech  tone  that  was  clearly  audible  ,  especially  for 
M>. 8.  A  resonating  field  between  the  nozzle  and  the  plate 
is  believed  to  be  responsible  for  that  screech  tone  re¬ 
ported  in  Refs.  7,  8  and  9.  The  resonance  phenomenon  is 
also  observed  in  the  normalized  power  spectra  of  Fig. 
(3-12).  The  energy  is  substantially  reduced  at  different 
frequency  bands,  and  is  Increased  in  a  narrow  band  around 
the  resonance  frequency. 

3.2.1  Surface  Pressure  Loading 

The  rms  of  surface  pressure  fluctuations  along  the 
plate,  normalized  by  the  dynamic  pressure  q  is  plotted  for 
M  ■  .5  and  .9  in  Fig.  (3-13).  At  Xg/d<5,  p'/q  has  a  peak 
at  rp/d»1.5.  This  is  the  approximate  location  where  vor¬ 
tex  rings  impinge  on  the  plate,  as  was  shown  earlier  in 
the  water  jet  (see  also  Ho  and  Kovasznay  (1974)].  For 
rg/d>3,  p'/q  is  independent  of  the  plate  location,  which 
is  a  characteristic  of  the  developed  wall  jet.  The  nor¬ 
malized  pressure  levels  arc  higher  for  M  ■  .9  th3n  for  M  - 
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.5  at  locations  rg/d<3.0.  This  is  shown  for  several  Mach 
numbers  in  Fig.  (3-14).  At  resonance,  M  *  .8  and  .9,  the 
surface  loading  is  generally  higher.  The  reason  for 
higher  pressure  levels  at  M  -  .8  than  at  H  «  .9  is  dis¬ 
cussed  later  in  connection  with  changing  the  plate  loca¬ 
tion  xQ/d. 

As  mentioned  before,  the  plate  is  divided  into  two 
regions,  the  impinging  region  and  the  wall  jet  region. 
Based  on  the  surface  pressure  fluctuations,  the  impinging 
region  is  further  divided  into  two  regions.'  an  inner  imp¬ 
inging  region  (0<rg/d<:.5)  and  an  outer  impinging  region 
(.5<rQ/d<2).  The  wall  jet  region  is  at  rQ/d>3.  These 
divisions  and  their  boundaries  are  rather  arbitrary  and 
are  made  for  convenience  in  the  future  discussion  of  dif¬ 
ferent  flow  characteristics.  The  effect  of  plate  location 
Xg/d  on  p'/q  for  these  regions  is  presented  in  Figs. 
(3-15)  and  (3-16).  In  the  inner  impinging  region,  p'/q 
increases  as  the  plate  moves  away  from  the  nozzle  (Fig. 
(3-15)].  Strong,  et  al.  (1967),  showed  that  p'/q  reaches 
its  maximum  value  of  0.12  at  Xg/d**7.  Fig.  (3-15)  also 
shows  that  p'/q  in  the  wall  region  is  independent  of  both 
M  and  Xg/d.  The  ^ow  there  is  so  dominated  by  the  wall 
through  the  viscous  forces  that  it  tends  to  "forget"  its 
past  history. 
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The  behavior  of  the  pressure  fluctuations  in  the 
outer  Inpinging  region  is  different  between  the  resonant 
and  nonresonant  jets.  In  the  resonant  case  of  M  -  .9, 
Fig.  (3-16a)  shows  that  the  pressure  fluctuations  p'/q 
goes  through  jumps  as  it  decreases  with  increasing  x^/d. 
These  jumps  are  similar  to  jumps  in  the  r "onance  frequen¬ 
cy  observed  when  x^/d  changes  as  reported  later  in  Sec. 
5.2  [also  see  Fig.  (5-4)].  At  M  «  .8,  Fig.  (3-16b) 
shows  that  p'/q  remains  at  a  value  .115  up  to  x^/d  -  2.7 
and  then  decreases  with  no  jumps  as  Xg/d  increases.  The 
reason  is  possibly  due  to  a  "weak"  resonance  at  M  ■  .8 
(see  Sec.  5.2).  The  surface  pressure  loading  in  the 
outer  impinging  region  is  due  to  the  impingement  of  the 
large  scale  structures  and  is  strongly  dependent  on  both 
the  jet  Mach  number  M,  and  the  plate  location  xQ.  The  va¬ 
lues  of  p'/q  seem  to  fall  between  two  limits  in  fig. 
(3— 16c):  a  lower  limit  of  a  non-resonant  case  at  M  «  .5 
and  an  upper  limit  of  the  resonant  case  at  M  *  .9.  The 
pressure  at  M  ■  .9  is  shown  to  jump  between  the  two  lim¬ 
its,  passing  through  an  intermediate  stage  with  values 
corresponding  to  the  "weak"  resonant  case  of  M  *  . 8.  The 
jumps  in  p'/q  at  M  *  .9  is,  therefore,  the  reason  for  the 
lower  levels  of  pressure  observed  in  Fig.  (3-14)  at  x^/d 
•  2.0  compared  to  the  M  «  .8  case. 
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3.2.2  Eddy-Convection  Velocity  on  the  Plate:  Comparison 
of  Different  Techniques 

The  convection  velocity  C  is  usually  evaluated  from 
the  correlation  function  R  by  either  of  two  ways  [Ro  and 
Kovasznay  (1976)]:  constant  separation  distance  ^  ^  with 
varying  the  time  delayT,  where  C  satifies  the  condition 

2R(  J,,  t) /  dr  &  o 

or  constant  time  delay  with  varying  the  separation  dis¬ 
tance  3  ,  where  C  satisfies  the  condition 

2R{  s  0 

The  calculated  values  by  the  two  methods  are  generally 
different.  Wills  (1964)  indicated  that  they  are  only 
equal  if  the  turbulent  flow  satisfies  Taylor's  hypothesis 
in  which  the  "frozen"  eddies  are  convected  unchanged.  To 
remove  the  ambiguity  in  evaluating  a  single  convection 
velocity,  he  argued  that  C  has  to  be  expressed  in  terms  of 
either  a  wave  number  k  or  frequency  f.  He  derived  the 
convection  velocity  C(k)  as  a  function  of  the  wave  number, 
and  defined  an  overall  convection  velocity  as  the  propaga¬ 
tion  velocity  of  the  peak  of  the  turbulent  energy,  integ¬ 
rated  over  different  wave  numbers  k.  This  method  is  suit- 


able  whenever  the  convection  velocity  of  different  wavel¬ 
ength  is  needed  as  in  water-wave  generation. 

In  the  present  study  the  frequency  is  chosen  as  a 

parameter  in  the  evaluation  of  convection  velocities.  The 

calculation  of  convection  velocity  at  any  frequency  C(f), 

especially  at  the  resonance  frequency,  is  more  appropriate 

in  studying  the  feedback  mechanism.  The  group  velocity 

C  (f)  is  defined  as  the  velocity  of  energy  propagation  at 
B 

frequency  f,  which  is  relevant  to  the  condition  used  by 
Wills  (1964)  to  evaluate  C(k).  Naturally,  for  nondisper- 
sive  waves  the  convection  velocity  C(f)  of  the  present 
study  and  C(k)  of  Wills  (1964)  are  the  same. 

In  this  section  the  convection  velocity  along  the 
plate  will  be  calculated  by  different  techniques  discussed 
earlier  in  Sec.  2.3.  Values  calculated  from  space-time 
correlation  functions  are  compared  with  those  using  con¬ 
stant-separation  correlation  functions.  The  narrow-band 
convection  velocity,  which  is  called  the  phase  velocity 
C(f),  evaluated  from  cro s s- s pe c tr a  function  is  also  pre¬ 
sented.  The  prewhitening  technique  is  used  to  calculate 
convection  velocities  for  resonant  cases. 

The  phase  velocity,  or  the  narrow-band  convection 
velocity  C(f)  is  calculated  from  the  relation 
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Cff )  *  <i ,  (3.1) 

where  ^(f)  Is  the  phase  angle  between  two  signals  at  po¬ 
ints  separated  by  The  group  velocity  is  proportional 

to  the  slope  of  the  phase  angle  and  is  given  by 

C,(f)  a  2  7Tj/ji£i£)  .  (5-2) 

3  / 

For  a  nondispersive  wave  the  slope  of  the  phase  angle  is 
constant  and  hence 

Cj  Of)  r  C(f)  s  constant  . 


An  example  of  the  cross-spectrum  and  phase  angle 
functions  between  two  points  on  the  plate  is  shown  in  Fig. 
(3  —  1  7a)  for  M  •  .8  and  xQ/d  -  6.  The  cross-spectrum  G(f) 
has  a  broad  peak  at  f*3200  Hz.  The  Strouhal  number  of 
this  peak,  fd/U  »  .30,  corresponds  to  the  frequency  of 
large  scale  structures  at  the  end  of  the  potential  core  in 
a  free  jet.  The  phase  angle  V(f)  fluctuates  around  a 
straight  line,  indicating  a  nondispersive  convected  wave. 
Misleading  data  appear  in  both  Cff)  and  -S^f)  st  low  fre¬ 
quencies  (f<700  Hz)  and  are  attr-.jted  to  the  low  frequen¬ 
cy  cutoff  in  evaluating  the  sp-c'-v”  *r  -  i  1. its  digital 
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record  length  [Eqn.  (2-1)].  Fig.  (3-1 7b)  shows  the 
coherence  function  Y(f)  for  the  seme  flow  parameters  of 
Fig.  (3-17a).  The  coherence  Tf(f)  takes  a  similar  shape 
as  G(f).  The  phase  angle  function  in  the  figure  is  evalu¬ 
ated  from  smoothed  frequency  components  as  mentioned  ear¬ 
lier  in  Sec.  2.3  (with  n  *  8).  The  phase  angle  function 
has  a  clear  constant  slope  over  the  shown  frequency  range. 
Eqn.  (3-1)  is  used  to  calculate  the  phase  velocity  from 
-?(f)  of  Fig.  (3-1  7b)  and  is  shown  in  Fig.  (3-1  7c).  As 
expected  from  the  constant  slope  of  -?(f),  the  phase  velo¬ 
city  is  the  same  for  different  frequencies  and  equal  to 
.62  of  the  Jet  exit  velocity  U,  which  is  indicative  of  the 
nond ispersiveness  of  the  convected  eddies  on  the  plate. 

At  low  jet  Mach  number  (M  •  .5),  Fig.  (3-18)  shows  a 
similar  cross-spectrum  function  and  a  constant-slope  phase 
angle  f  unction  . 

The  correlation  coefficients  are  presented  in  Figs. 
(3-19)  and  (3-20).  The  correlation  coefficients  are  plot¬ 
ted  as  a  function  of  separation  distance  and  time  delay 
T for  the  case  of  M  ■  . 8  in  Fig.  (3-19).  The  characteris¬ 
tic  forms  of  these  correlation  curves  are  in  agreement 
with  those  of  Refs.  3  and  6.  The  correlation  coefficient 
wa s  very  small  (<5X)  between  any  point  in  the  inner  imp¬ 
inging  region  and  the  rest  of  the  plate.  This  is  possibly 
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due  to  turbulent  flow  circulation  in  the  "atagnation  bub¬ 
ble"  reported  by  Donalson,  et  al »  (1971)*  Point  2  (r^/d 
-  1.0)  ie  choeen  as  a  fixed  reference  point.  This  point 
is  in  the  region  where  the  vortex  rings  were  observed  to 
impinge  on  the  plate  in  the  water  jet  experiment.  The 
figure  shows  a  50X  correlation  coefficient  which  is  almost 
constant  over  the  outer  impinging  region  (l<rQ/d<2).  The 
correlation  coefficient  drops  for  larger  separation  dis¬ 
tance.  since  convected  eddies  tend  to  lose  their  identi¬ 
ties  due  to  viscous  dissipation  and  the  development  of 
small  scale  turbulence. 

At  low  jet  speed  (M  -  .5),  the  correlation  coeffi¬ 
cient  was  found  to  drop  sharply  for  a  separation  distance 
larger  than  one  nozzle  diameter.  The  reason  is  that  large 
scale  structures  are  less  coherent  for  a  weak  feedback  at 
such  low  jet  speed.  The  correlation  function  is,  there¬ 
fore,  shown  in  Fig.  (3-20)  as  a  function  ofTwith  con¬ 
stant  separation  distance  equal  to  d/2.  The  correlation 
coefficient  between  two  points  in  the  inner  impinging  re¬ 
gion  Rn  (T)  is  again  small  and  of  the  same  order  of  mag- 
nitude  as  the  noise.  It  increases  as  we  enter  the  outer 
impinging  region  where  the  large  coherent  structures  im¬ 
pinge  on  the  plate. 

At  high  subsonic  jet  speed  and  small  nozzle  to  plate 


distance,  resonance  occurred  and  the  pressure  signal  was  a 
distorted  sine  wave.  The  correlation  function  also  took  a  ; 
sinusoidal  pattern  with  no  clear  optimum  peak  to  evaluate 
eddy  convection  velocity.  The  prewhitening  technique  re¬ 
ported  in  Sec.  2.3  was  used. 

1  e  raw  surface  pressure  signals  at  different  loca¬ 
tions  on  the  plate  for  a  resonant  Jet  at  M  ■  .9  and  Xg/d  * 
4.5  is  shown  in  Fig.  (3-21).  The  signal  at  rg/d  -  0.  is 
almost  a  pure  sine  wave  owing  to  the  induced  field  of  the 
Impinging  coherent  eddies.  This  is  supported  by  the  rela¬ 
tively  low  rms  value  of  the  pressure  fluctuations  at  r^/ d 
■  0.  in  Fig.  (3-13).  The  signals  at  the  outer  impinging 
region  exhibit  rather  distorted  sine  waves  because  of  the 
small  scale  turbulence. 

The  power  spectra  of  these  signals  are  presented  in 
Fig.  (3-22).  Most  of  the  energy  is  concentrated  in  a 
narrow  band  around  the  resonance  frequency  at  the  location 
rg/d  *  0.  In  the  outer  impinging  region  other  frequency 

components  appear  in  the  spectra.  It  is  noticeable  that 
at  r q / d  *  .5,  the  energy  is  more  evenly  distributed  over  a 
wide  band  of  frequencies  due  to  the  random  flow  in  the 
"stagnation  bubble".  In  all  the  spectra,  a  smaller  peak 
is  observed  at  the  first  harmonic  of  the  resonance  fre¬ 
quency. 
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The  cross-spec trum  between  two  points  is  shown  In  ' 

| 

Fig*  (3-23a).  Two  peaks  at  the  resonance  frequency  and  j 
its  first  harmonic  is  shown*  A  narrow-band  digital  filter  I 
(bandwidth  «  312  Hz)  was  used  in  the  prewhitening  process 
to  average  out  these  two  peaks*  The  resr’ t  is  shown  in 
Fig*  ( 3  —  2 3 b )  with  the  corresponding  phase  angle  function* 
The  real  gain  from  prewhitening  is  shown  in  the  correla¬ 
tion  functions  of  Fig.  (3-24).  As  expected,  the  correla¬ 
tion  functions  before  prewhitening  are  sinusoidal  with 
respect  to  the  delay  time  ,  while  the  prewhitened  corre¬ 
lation  function  shows  a  distinctive  peak.  The  delay  time 
of  this  peak  is  used  to  evaluate  a  broad-band  convection 

velocity  C  for  resonant  cases.  Fig.  (3-25)  shows  a 

space-time  prewhitened  cross-correlation  functions  of  a 
resonant  case  (M  ■  .9  and  x^/d  “  4.5).  The 

cross-correlation  function  for  the  inner  impinging  region 
is  plotted  separately  in  Fig.  (3-25a).  The  correlation 
between  point  0  and  1  shows1 a  negative  time  delay  peak, 
indicating  a  reversed  flow  towards  the  stagnation  point. 
This  is  consistent  with  flow  circulation  in  the  stagnation 
bubble  as  reported  by  Donalson,  et  al.  (1971).  The 

correlation  coefficients  in  the  outer  impinging  region  are 
shown  in  Fig.  (3-25b).  Here,  the  correlation  coeffi¬ 
cients  show  a  continuous  decrease  with  Increasing  the  sep¬ 
aration  distance  and  time  delay,  which  is  different  from 
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.  the  behaviour  of  those  in  Fig.  (3-19).  The  reason  is 
!  that  the  coherent  part  of  the  signals  used  to  produce  the 
correlation  functions  has  been  reduced  substantially  dur¬ 
ing  the  prewhltening  process. 

The  convection  velocity  C  along  the  plate  is  plotted 
in  Fig.  (3-26).  Except  for  a  negative  value  near  the 
stagnation  point,  the  convection  velocity  Increases  to  a 
peak  at  r^/ d  ■  1.25  and  then  decreases  asymptotically  to  a 
value  of  .45  U.  The  higher  convection  velocity  near  the 
outer  impinging  region  (*.65  U)  agrees  with  observation  of 
pictures  and  movies  [see  Figs.  (3-3)  and  (3-7)]  of  the 
impinging  water  jet.  The  shear  layer  was  observed  to 
curve  as  it  approaches  the  plate.  The  deceleration  is 
balanced  by  the  acceleration  due  to  shear  layer  curvature, 
and  the  vortex  rings  were  seen  to  maintain  their  free 
speed  (or  even  accelerate)  during  impingement. 

The  data  points  in  Fig.  (3-26)  for  M  ■  .  8  are  evalu¬ 
ated  from  the  separation  distance  and  the  optimum  time 
delay  of  two  different  correlation  plots.  The  convection 
velocity  C  is  plotted  at  the  half  point  of  the  separation 
distance.  The  values  of  C  at  r^/d  »  1.75  and  2.5  are  eva¬ 
luated  from  a  constant  separation  distance  correlation 
function  (withlj"  d/2)  similar  to  those  of  Fig.  (3-20). 
These  points  are  in  accord  with  others  calculated  by  the 


variable  separation  distance  correlation  functions  of  Fig* 
(3-19).  As  previously  noted,  this  should  not  happen  un¬ 
less  Taylor's  hypothesis  is  satisfied.  Lin  (1952)  has 
shown  that  Taylor's  hypothesis  is  valid  only  if  the  turbu¬ 
lence  level  is  low,  viscous  forces  are  negligible,  and  the 
mean  shear  is  small.  It  is,  therefore,  unreasonable  to 
expect  the  hypothesis  to  apply  to  the  flow  of  an  impinging 
Jet  on  a  plate.  Hence,  the  concept  of  convected  "frozen" 
eddies  is  excluded  as  a  reason  for  the  agreement  in  the 
velocity  calculated  by  the  two  methods.  To  understand 
this  result,  the  expected  difference  in  the  calculated  va¬ 
lues  of  C  by  the  two  correlation  function  plots  must  be 
discussed  first.  As  the  separation  distance  Increases, 
small  eddies  are  bound  to  lose  their  identity  and  become 
uncorrelated  in  comparison  with  larger  size  eddies.  This 
is  in  effect  a  lowpass  filter  with  a  cutoff  frequency  de¬ 
creasing  with  increasing  added  during  the  evaluation  of 
correlation  coefficients.  With  a  fixed  the  correlation 
estimates  are  most  sensitive  to  the  convection  of  eddies 
with  a  size  corresponding  to  the  fixed  bandwidth  filter. 
The  consistency  of  the  calculated  convection  velocity  ir¬ 
respective  of  the  effective  filter  used  indicates  that  ed¬ 
dies  with  different  sizes  convect  at  the  same  speed.  In 
other  words,  if  one  uses  the  method  proposed  by  Wills 
(1964),  the  calculated  convection  velocity  for  different 
wave  numbers  C(k)  would  be  constant.  Thjs,  the  flow  over 


It  was  shown  in  the  pictures  of  Sec.  3.1  that  the 
coherent  structures  have  a  vortex  ring  shape.  At  high  Re¬ 
ynolds  number  (Re>8000)  a  higher  instability  mode,  namely 

I 

a  helical  mode,  appeared.  Browand  and  Laufer  (1975)  ob- 

f 

served  both  the  axisymmetric  and  helical  modes  in  a  free 

Jet  at  Re  ■  5000  to  20,000.  Fuchs  and  Michel  (1977)  used  ■ 

i 

a  set  of  near  field  microphones  placed  in  a  plane  perpen-  j 

* 

dicular  to  and  at  equal  distances  from  the  jet  axis.  By  ! 

• 

Fourier  analyzing  the  real  part  of  the  cross-spectrum 

between  these  microphones,  they  were  able  to  determine  the 

azimuthal  constituents  of  the  turbulent  energy  at  any  fre- 

4  7 

quency.  For  a  range  of  10  <Re<2xl0  (.3  <  M  <  .7),  they 
found  that  the  energy  was  mostly  contained  in  the  axisym¬ 
metric  mode,  although  higher  modes  were  present.  Neuvorth 

(1973)  observed  that,  while  at  M  ■  . 5  pictures  showed  an 
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axisymmetric  coherent  structures,  at  M  ■  .9  a  helical  mode 
appeared  In  the  free  jet.  However,  when  the  plate  was  In¬ 
serted  his  pictures  and  those  of  Vagner  (1971)  showed  that 
only  the  axisymmetr ic  mode  remained.  Neuworth  predicted 
that  waves  reflecting  from  the  plate  were  responsible  for 
the  axisymmetric  excitation  of  the  free  shear  layer  near 
the  nozzle. 

In  the  present  study  the  stability  modes  were  exam¬ 
ined  by  placing  five  pressure  transducers  at  equal  dis¬ 
tance  from  the  stagnation  point  on  the  plate.  Gross  fea¬ 
tures  of  the  stability  modes  of  the  impinging  jet  were  ex¬ 
amined  from  the  correlation  function  between  points  on  the 
plate  at  different  azimuthal  angles.  An  example  of  the 
correlation  function  between  two  points  separated  by  * 
270  is  shown  in  Fig.  (3-27).  The  correlation  for  this 
resonant  case  has  a  maximum  at T ■  o,  which  suggests  that 
the  axisymmetric  mode  is  dominant.  The  prewhitened  corre¬ 
lation  function  is  also  shown  with  a  peak  atT*  0,  indi¬ 
cating  that  the  turbulent  portion  of  the  surface  pressure 
fluctuations  are  also  dominated  by  an  axisymmetric  mode. 
Similar  results  are  shown  for  different  azimuthal  angles 
in  Fig.  (3-28).  Furthermore  Fig.  (3-29)  shows  that  the 
phase  angle  is  approximately  equal  to  zero  over  all  fre¬ 
quencies,  which  suggests  that  the  existing  turbulent 
scales  are  axisymmetric. 

55 


Similar  analysis  was  carried  out  to  ensure  axisymme- 
try  in  the  flow  field  between  the  nozzle  and  the  plate* 
Fig.  (3-30)  shows  the  cross  correlation  function  before 
and  after  prewhitening  between  two  transducers  at  the  exit 
plane  of  the  nozzle.  Axisymmetry  is  indicated  by  the 
peaks  at  T «  0. 
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CHAPTER  4 

THE  TWO  BRANCHES  OF  THE  FEEDBACK  LOOP 


In  the  previous  chapter,  general  characteristics  of 
the  impinging  jet  were  presented.  The  validity  of  differ¬ 
ent  digital  data  techniques  was  examined.  The  convection 
velocity  on  the  plate  was  found  to  be  consistent  when  eva¬ 
luated  by  these  techniques,  and  to  agree  well  with  results 
from  Refs.  3  and  6  .  Finally,  the  flow  field  of  the  imp¬ 
inging  jet  was  found  to  be  dominantly  axisynmetric  . 

As  mentioned  earlier  in  Chapter  1,  a  conjecture  for 
the  feedback  mechanism  causing  the  resonance  could  take 
the  following  form:  coherent  eddies  are  convected  down¬ 
stream  and  impinge  on  the  plate;  waves  are  generated 
which  propagate  upstream  and  excite  the  shear  layer  near 
the  nozzle;  the  shear  layer  rolls  into  new  coherent  ed¬ 
dies  and  the  feedback  loop  is  closed.  In  this  chapter, 
pressure  measurements  are  used  to  detect  waves  propagating 
in  the  near  field  between  the  notsle  and  the  plate.  The 
characteristics  of  tv v—vtf  are  studied  and  the  rela¬ 
tion  between  them  at  :  i  i:'  i  rve  r : ;  ca  ted  . 
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4.1  CORRELATIONS  OF  NEAR  FIELD  PRESSURE  MEASUREMENTS 


Power  spectra  of  signals  from  near  field  microphones 
are  shown  in  Fig.  (4-1)  for  M  «  .  8  and  Xg/d  ■  7.  The 
spectra  do  not  show  a  distinguished  high  peak,  indicating 
no  resonance.  The  correlations  of  these  signals  are  plot¬ 
ted  in  Fig.  (4-2)  as  functions  of  separation  distance  and 
time  delay  T.  For  positiveT,  the  correlation  coefficient 
peaks  indicate  a  downstream  convected  wave.  On  the  other 
hand,  other  peaks  at  negative  time  delays  are  also  pre¬ 
sent.  While  the  peak  in  R  (T)  at  negative  time  delay 
may  appear  as  a  result  of  normal  oscillations  in  the 
correlation  function,  the  one  in  R  (T)  does  not.  These 
peaks,  therefore,  indicate  upstream  propagation,  and  the 
behavior  of  the  correlations  suggest  the  possibility  of 
two  waves  propagating  in  the  near  field. 


At  resonance  the  correlation  functions  are  sinusoidal 
with  respect  to  T  and  provide  no  information  about  wave 
propagation.  However,  the  presence  of  the  two  waves  is 


observed  in  the  cross-spectra  and  phase  angle  functions  in 
Fig.  (4-3).  The  cross-spectrum  between  two  microphones 
near  the  axis  of  the  jet  G  (f)  is  dominated  by  two 

i#  u 

strong  resonance  peaks.  The  presence  of  two  resonance 


peaks  instead  of  one  is  attributed  to  the  tendency  of  the 
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jet  to  intermittently  jump  between  two  resonance  frequency 
stages  in  this  specific  case  (M  ■  .9  and  x^/d  ■  5)*  These 

t 

frequency  Jumps  are  discussed  later  in  Sec.  5.2.  The  ' 

1 

phase  angle  ^  (f),  despite  some  discrepancies  at  5.5  1 

I»3I 

KHz<  f  <6  KHz,  increases  with  frequency.  These  discrepan¬ 


cies  may  be  due  to  very  low  levels  of  energy  at  this  fre- 

-4 

quency  band  [G  ( f  >  < 1 0  ].  Here  the  noise  signal  is  ex- 

**n 

pected  to  be  of  the  sane  order  of  magnitude  as  the  meas¬ 


ured  signals.  However,  the  positive  slope  of  (f)  in- 
dlcates  a  downstream  convected  wave.  Although  G  (f) 
shows  peaks  at  the  same  resonance  frequencies,^  (f)  be- 
haves  differently.  The  negative  slope  of  ^  (f)  for 
these  microphones  located  further  away  from  the  jet  axis 


indicates  the  existence  of  an  upstream  wave. 


The  cross-correlations  between  microphones  I  &  33  be¬ 
fore  and  after  prewhitening  are  shown  in  Fig.  (4-4)  and 
the  corresponding  cross- cor  relations  between  HI  &  S'  are 
shown  in  Fig.  (4-5).  The  correlation  function  R_  (T)  is 
a  sine  wave  as  expected  from  two  resonant  signals. 
However,  the  prewhitened  correlation  R  (T)  shows  two  op- 
timum  peaks  at  positive  and  negative  time  delays.  This  is 
in  agreement  with  earlier  observations  in  the  non-resonant 
case  of  Fig.  (4-2).  In  Fig.  (4-5)  the  correlations 

V 

R_  (T)  and  R__  (TO  measured  at  locations  further  away 

Ju»ar  nr,  it 

from  the  jet  axis  show  only  negative  tic  >.  delay  peaks,  in- 


dicating  predominant  upatream  propagation. 

The  above  two-point  measurements  provide  evidence  for 
the  presence  of  two  waves  propagating  in  the  near  field. 
Furthermore,  the  upstream  propagating  wave  is  shown  to  be 
dominant  in  the  near  field  away  from  the  jet  axis.  This 
provides  an  opportunity  to  study  one  of  the  two  waves, 
namely  the  upstream  wave,  through  measurements  taken  away 
from  the  jet  axis. 

4.1.1  Characteristics  of  the  Upstream  Wave 

The  pressure  fluctuations,  p/  ,  measured  in  decibels 
by  a  microphone  moving  outwards  in  a  plane  perpendicular 
to  the  jet  axis  at  the  nozzle  exit  is  plotted  in  Fig. 
(4-6).  The  figure  shows  p  “  for  a  free  jet  at  H  »  .8  com¬ 
pared  to  an  impinging  jet  at  M  -  .8  and  .9.  In  a  free 

«\ 

jet,  p'  drops  sharply  from  r/d  -  .5  to  r/d  ■  1.0  (not 
shown  in  the  figure)  and  then  descends  at  a  rate  propor¬ 
tional  to  r”*’^.  This  is  in  agreement  with  Chu  (1975), 

who  also  showed  that  the  near  field  of  the  present  jet  ex- 

o  2 

tends  to  r/d-  30  at  ©  -  90  .  For  the  impinging  jet  p 

has  higher  levels,  but  behaves  like  the  free  jet  up  to  r/d 
-  3.  These  high  levels  are  mainly  due  to  added  energy 

from  the  upstream  wave.  At  large  radii,  p'  increases  for 
3  <r/d<  6.5  and  the.',  decreases  at  a  rate  proportional  to 


_2 

r  .  This  indicates  that  the  energy  is  doninated  by  the 
upstream  wave  for  r/d>3,  and  that  the  far  field  of  this 
wave  is  reached  at  r  *  6.5,  after  which  p*  follows  the 
Inverse  square  law. 

The  phase  difference  at  the  resonance  frequency 

between  two  microphones  in  a  plane  perpendicular  to  the 

jet  axis  at  the  nozzle  exit  is  shown  in  Fig.  (4-7).  One 

of  the  microphones  was  fixed  near  the  nozzle  lip  at  r/d« 

.5,  while  the  other  was  allowed  to  move  radially  away  from 

the  jet  axis.  For  r/d>3,  the  phase  angle  increases  almost 

linearly  with  r.  In  this  region  the  upstream  wave  was 

shown  to  be  dominant  from  the  pressure  measurements  of 

Fig.  (4-6),  and  also  from  the  spectra  and  correlation 

functions  of  Figs.  (4-3)  and  (4-5).  Therefore,  the  phase 

angle  plotted  against  r  for  r/d>3.0  represents  an  almost 

plane  wave  front  propagating  in  the  upstream  direction. 

The  direction  of  propagation  is  the  perpendicular  to  the 

average  slope  d^  (f  ) /dr 

r 

For  r/d<3.0,  Fig.  (4-7)  shows  that  the  two  signals 
are  in  phase  except  for  a  "bulge"  around  r/d  ■  2.  This 
bulge  is  believed  to  be  an  erroneous  deviation  from  a  zero 
phase  value  (a  zero  phase  at  the  nozzle  exit  is  discussed 


ponents  [Fig*  (4-8)],  and  that  the  phaae  angle  wa a  calcu¬ 
lated  at  a  single  frequency  f  by  averaging  over  a  long 
tine  record*  The  appearance  of  random  high  frequencies  in 
the  signals  at  r/d  -  2  is  possibly  due  to  the  following 
two  reasons: 

1.  The  acoustic  waves  radiated  from  the  stagnation 
region  of  the  plate  to  this  region  suffer  (as  will  be 
shown  later)  from  random  scattering  as  they  travel  through 
a  thick  turbulent  shear  layer,  and/or 

2.  A  microphone  in  this  vicinity  is  sensitive  to 
pressure  integrated  over  the  whole  flow  field. 

Fig.  (4-8)  shows  that  a  resonant  sinusoidal  signal 
dominates  close  to  the  edge  of  the  shear  layer  at  r/d«*.5, 
and  away  from  the  jet  axis  at  r/d  ■  6.5  where  the  upstream 
wave  is  strong 

* 

In  order  to  evaluate  the  speed  of  the  upstream  wave, 
a  third  microphone  JJJ  is  placed  downstream  of  the  moving 
one.  The  phase  angled  (f)  is  plotted  in  Fig.  (4-9)  as 

m>ix 

a  function  of  frequency  f.  Three  samples  of  (f)  at 
different  / d  are  shown  in  the  figure.  The  phase  angle 
increases  linearly  with  f.  Since  M  is  delayed,  the  con¬ 
stant  positive  slope  i'  ’.'.-tes  a  nond  ispersive  upstream 


*  a 


propagating  wave.  It  is  also  noticeable  that  fluctuations 
in  0  (f)  around  Its  mean  decrease  as  r_,  /d  Increases. 

nr, a  n 

owing  to  the  dominant  contributions  of  the  upstream  wave 
over  the  pressure  signals  in  the  near  field  away  from  the 
jet  axl6.  It  is  now  clear  that  the  upstream  wave  can  be 
treated  as  a  nondispersive  [Fig.  (4-9)],  plane  [Fig. 
(4-7)]  wave.  Furthermore,  the  angle  Q  between  the  jet 
axis  and  the  direction  of  propagation,  can  be  calculated 
from  the  mean  slope  of  the  wave  front  according  to  the  re¬ 
lation 


=  w 


-I 


<Mfr) 


(4-1) 


From  Fig.  (4-7)  a  value  of©a*3$  is  calculated  using 
Eqn .  (4-1).  This  value  is  further  confirmed  from  the 

prevhitened  cross  correlation  functions  between  II  and  HT . 
The  correlation  coefficients  optimum  peaks  R_  (T* )  are 

m,n  p 

plotted  versus  r^.  /d  i.n  Fig.  (4-10).  A  broad  peak  of 
R  (7* )  shown  around  r_  /d  «  3.3  corresponds  to  a  propa- 

nr.tr  p 

gati  n  direction  33  .  This  is'  in  good  agreement  with 


the  value  calculated  above. 


The  speed  of  the  nondispersive  upstream  wave  C.,  can 

* 

be  calculated  from  the  slope  of  the  phase  angle  between 
two  microphones  (e.g.  31  &JU  )  by  the  following  relation 
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where  ©a  is  the  unit  vector  in  the  direction  of  upstream 
wave  propagation  inclined  by  an  angle  ©a  to  the  jet  axis  . 
Eqn •  (4-2)  is  used  to  plot  the  magnitude  of  the  upstream 

wave  speed  normalized  with  respect  to  the  ambient  speed  of 
sound  for  different  r^  / d  in  Fig.  (4-11).  It  is 

shown  that  the  upstream  wave  travels  with  a  speed  that  is 
equal  to  the  speed  of  sound  a^. 

» 

A  constant  value  of  © 35  from  the  mean  slope 
dy(fr)/d(r/d)  of  Fig.  (4-7)  was  used  in  Fig.  (4-11). 
However,  Fig.  (4-7)  shows  that  the  wave  front  is  slightly 
'  curved  and  that  d«?(  1  )  /  d  (  r  /  d)  Increases  as  r/d  increases. 
This  resulted  in  smaller  values  of  C^/&q  for  larger  r^.  /d, 
which  is  the  trend  one  observes  in  the  figure. 

The  direction  of  the  upstream  wave  propagation  for 
different  plate  locations  is  presented  in  Fig.  (4-12). 
It  is  shown  that  as  the  plate  moves  closer  to  the  nozzle, 
the  upstream  wave  propagates  at  a  larger  angle  relative  to 
the  Jet  axis. 
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4.2  A  SIMPLE  MATHEMATICAL  MODEL 

The  correlation  and  cross-spectra  functions  of  the 
measured  pressure  signals  have  shown  a  downstream  and  an 
upstream  propagating  wave.  The  upstream  wave  can  be  tre¬ 
ated  as  a  plane  sound  wave  propagating  along  a  path  in¬ 
clined  at  an  angle  0a  to  the  jet  axis.  The  following  sim¬ 
ple  mathematical  model  presupposes  two  coun terpropagating 
waves  in  the  near  field  of  the  Jet.  This  model  will  help 
to  explain  the  feedback  mechanism  through  analysis  of 
measured  correlation  functions. 


The  pressure  at  any  point  in  the  near  field  between 
the  nozzle  and  the  plate  is  the  superposition  of  two  parts 


f>(tj  f  )  =  oC  (t  )  f  (tj  }  )  j  (4-3) 


where , 

5  is  the  position  vector  with  respect  to  the  center 
of  the  nozzle  exit  plane  (see  Fig.  (4-13)],  and  o<  and  are 
the  two  parts  of  the  signal  due  to  the  downstream  and  the 
upstream  travelling  waves,  respectively.  Eon.  (4-3)  is 
also  expressed  as  a  pair  of  travelling  waves,  i.e. 


pa,})  =ali)e  + 


& \ S)  €  j 
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where  •<J>.  b(J)  are  the  wave  amplitudes, 

kj,k^  ere  the  downstream  and  the  upstream  wave  numbers 
respectively,  and 

Is  a  constant  phase  shift  between  the  two  waves. 


The  following  assumptions  are  made  In  order  to  sim- 
plify  the  model: 


1.  The  two  waves  are  assumed  to  be  plane  waves.  It 

Is  anticipated  that  the  downstream  wave,  which  is  due  to 

the  convected  large  scale  structures,  Is  travelling  along 

the  jet  axis.  The  upstream  wave  propagates  In  a  direction 

making  an  angle  0  with  the  jet  axis. 

ft 


2.  The  waves  are  monochromatic,  i.e.  only  one  fre¬ 
quency  component  at  the  resonance  frequency  (or  the  domi¬ 
nant  frequency  in  a  non-resonant  case)  is  considered.  The 
two  waves  are  also  assumed  to  have  the  same  frequency 
since  the  upstream  one  is  generated  from  the  incidence  of 
the  downstream  wave  on  a  solid  surface,  i.e. 


uj,(k()  -  Vj{kz)  ■=  u>r 


Eqn .  (4-4)  then  takes  the  fo~r 


P  = 


.-Vt-f 


H-5) 


By  definition,  the  phase  velocity  of  the  downstream  wave 
in  the  x-direction  is 


(4-4) 


and  the  phase  velocity  of  the  upstream  wave  is 


C2  -  _j£r.  €L 

IkJ  ^ 


(4-  7) 


The  phase  shift  in  Eqn.  (4-5)  can  be  evaluated  by 
applying  the  boundary  condition  at  the  plate.  Since  the 
upstream  wave  is  assumed  to  be  generated  by  the  downstream 
wave  as  it  impinges  on  the  plate,  the  phase  difference 
between  the  two  waves  should  be  equal  to  zero  there.  The 
phase  difference  between  the  two  waves  at  any  location  is 


9  =  -  k,  x  _  k2.  f  +  *9q  . 


(4-8) 


The  boundary  condition  is 


VcO  c it  1  =  J 


S  • 


where  is  the  position  vector  of  an  apparent  sound 
source  on  the  plate  fFi;-..  (4-13)].  Substituting  (4-9)  in 

(4-8 ) ,  one  obtains 


*0  s  *0  •*  Ki  •  • 


( 4-/o  ) 


Hence  eqn.  (A-8)  may  be  written 


*  =  K,(x0-x)  +  K 2*ns-1)  >  H-H ) 


which  represents  the  phase  difference  between  the  two 
waves  at  any  location  in  the  near  field  between  the  nozzle 
and  the  plate* 


A. 2.1  The  Correlation  Function 


From  Eqn.  ( A — 5 )  the  pressure  at  points  I  and  Hare 


i(«4-K.*,4-?0)  cKt  +  Vft) 

Pj  =  Q| e  •+  b,  e 


=  Ql  «, 


b,  P, 


i(«A-MlrtS’0)  L^rt+kj.f  ) 

p  =  v  +  b*  «• 


(4-/2) 


-  Q2 


^2  5 


O 
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where  oi  j  ,  (&  j,  oL  2*  ft  2  are  the  si8nal  conponents  due  to  the 
two  waves. 

The  correlation  function  of  two  signals,  each  is  the 
sun  of  two  stationary  processes  is  given  in  Ref  30. 
Similarly,  the  correlation  function  of  the  two  signals  of 
Eqn.  (4-12)  is 


(4-13) 


The  above  equation  shows  the  correlation  function  (70 
as  the  algebraic  sura  of  the  individual  component  correla¬ 
tion  functions.  The  correlation  functions  of  these  compo¬ 
nents  are  by  definition: 


R  (T)  =  ±  e 


«K r-MXn-*!)] 


W'i  e 


i[o,T+  j 

A 


t£o>rT“  k,  xn  “  ^2’^+  ^0] 
R  (TJ  -±C 
*0*  2 


R  (TJ  rj_ 

•ift  * 


(4-Ha) 

(-<-.14  b) 

(4  -MO 

(4-  I4d) 


Using  Eqn.  (4-6),  Eqn.  (4-1 4a)  gives  an  optimum  tim< 
delay 


(I) 

T  * 


Xw  —  Xi 


(4-  IS) 


and  using  Eqn.  (4-7),  Eqn.  (4—1 4b)  similarly  gives 


-  Vn  •  4 


(4-16) 


where  Cj  is  the  speed  of  the  upstream  wave  inclined  at  ar 
angle  ©  to  the  jet  axis.  Similarly  from  Eqns.  (4-10), 

a 

(4-1 4 c)  ,  and  (4-14d) 


■-[ 


*0  ~  *1T  +  %T,S  * 

c, 


(-♦-17) 


4 


Xp  -  X 


kg-  e, 


i-  L  - 1«) 


where  is  the  vector  connecting  I  and  S,  and  Y  is 

M  8  ^SpII 

the  vector  connecting  S  and  II  [Fig.  (4—13)5* 


Eqns.  (4-13)  and  (4-15)  to  (4-18)  predict  that  the 
correlation  function  between  two  points  in  the  near  field 


'Vf* 


exhibits  two  peaks  with  positive  time  delays  as  well  as 
two  peaks  with  negative  time  delays.  The  physical  in¬ 
terpretation  of  the  four  time  delays  of  these  peaks  is  ex-  i 

plained  as  follows:  t  is  the  time  it  takes  the  down- 

u) 

stream  wave  to  travel  from  I  to  II;  X  is  the  time  it 

takes  the  upstream  wave  to  travel  from  II  to  I  and  should 

(4) 

appear  as  a  negative  tine  delay  in  R  (f)  ;  7“  is  the 

IjU 

time  it  takes  the  downstream  wave  to  travel  from  I  to  the 

plate  with  velocity  Cj,  plus  the  time  it  takes  for  the 

generated  wave  to  propagate  back  from  the  plate  to  II; 
(3) 

T  is  the  time  for  the  downstream  wave  to  go  from  II  to 

the  plate  plus  the  time  for  the  generated  upstream  wave  to 

<C3) 

bounce  back  to  I.  C  should,  therefore,  appear  as  a  ne¬ 
gative  time  delay  in  R  (T) . 


The  correlation  function  between  the  near  field  pres¬ 
sure  signal  and  a  surface  pressure  signal  in  the  impinging 
region  on  the  plate  can  be  derived  by  substituting 

x  *xA  and  X  ■ f  in  Eqns.(4-15)  to  (4-18)  which  yields 

H  u  H  s 


<•>  (4)  *  * 

T  =  T  r  *0  z  *  1  > 

c, 


<41  (3) 

T  s  T  - 


c. 


La.xq) 
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These  leter  results  show  that  the  four  peaks  degenerate 
Into  only  two  In  the  correlation  between  a  point  I  in  the 
near  field  and  a  point  S  on  the  plate.  The  positive  time 
delay  of  one  peak  (Eqn.  (4-19)]  corresponds  to  the  down¬ 
stream  convection  to  the  plate  and  the  negative  time  delay 
of  the  other  peak  tEqn.  (4-20)]  corresponds  to  the  up¬ 
stream  propagation  from  the  plate  to  point  I. 


For  the  special  case  of  x^»  x^-  x^,  Eqn.  (4-13) 

gives  the  autocorrelation  function  R.  .  (*r)  with  optimum 
time  delays  from  Eqns.  (4-15)  to  (4-18)  given  by 


(O  U ) 

T  r  s  0 


((4-20 


<*) 

r-  _r  = 


where  the  index  i  refers  to  positions  I,  II,  III,  or  IT  . 
Eqn.  (4-21)  predicts  that,  due  to  the  presence  of  the 
plate,  an  autocorrelation  function  should  exhibit  an  extra 
peak  besides  the  usual  one  atT*  0.  The  time  delay  of 
that  extra  peak  corresponds  to  the  time  it  takes  the  down¬ 
stream  wave  to  travel  to  the  plate  combined  with  the  time 
it  takes  for  the  generated  vave  to  propagate  back  to  the 
same  point.  Eqn.  (4-21)  £.’=■'  indicates  that  the  extra 
peak  appears  in  both  the  e  and  positive  time  dela- 
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ys,  which  is  consistent  with  the  autocorrelation  being  an 
even  function. 

4.3  MEASURED  CORRELATIONS:  COMPARISON  WITH  THE  MODEL 

The  pressure  in  Eqn .  (4-5)  is  expressed  as  the  sum 

of  two  waves  that  have  the  same  single  frequency  component 
u>r«  The  frequency  is  the  dominant  resonance  frequency 
or  the  peak  frequency  of  the  narrow-band  signals  for 
non-resonant  cases.  However,  the  measured  pressure  sig¬ 
nals  contain  other  frequency  components  due  to  the  random¬ 
ness  of  the  turbulent  flow.  This  randomness  adds  unique 
characteristics  to  each  convected  event  in  the  flow  field. 
Accordingly,  the  predicted  optimum  time  delays  of  Eqns. 
(4-15)  to  (4-18)  are  expected  to  correspond  to  peaks  of 
high  correlation  levels  in  the  correlation  coefficients  of 
measured  signals.  The  correlation  functions  for  the 
non-resonant  cases  are  shown  in  Figs.  (4-14)  to  (4-16), 
and  the  prewhitened  correlation  functions  for  a  strong  re¬ 
sonant  case  are  displayed  in  Fig.  (4-17). 

The  autocorrelation  functions  in  Fig.  (4-14)  for  M  * 
.8  and  / d  ■  7  show  extra  peaks  [peaks  (3)  and  (4)]  be¬ 

sides  the  usual  unity  peak  atT  ■  0,  in  agreement  with  the 
model.  The  space-time  correlations  for  this  non-resonant 

case  are  plotted  in  Fig.  (4-15).  For  zero  separation, 
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R_  _ (t)  shows  that  T  ■  T  ■  0.  As  the  separation  dis- 
i » I 

tance  increases,  peak  (1)  moves  to  larger  positive  time 
delays  Indicating  a  downstream  convection,  while  peak  (2) 
moves  in  the  opposite  direction  indicating  an  upstream 
convection.  Peaks  (3)  and  (4)  are  at  larger  time  delays 
and  they  meet  with  (1)  and  (2)  in  the  correlation  function 
with  the  transducer  on  the  plate  R_  n(f)  ,  as  predicted  by 
Eqns.  (4-19)  and  (4-20). 

The  approximate  loci  of  these  peaks  form  the  letter 
"W"  as  indicated  by  the  dashed  lines  in  Fig.  (4-16).  In 
this  figure,  the  correlation  functions  are  separated  vert¬ 
ically  in  proportion  to  actual  physical  separation  dis¬ 
tance  between  microphones,  along  the  jet  axis.  For  a  con¬ 
stant  convection  velocity  in  the. x-d irec t ion ,  the  dashed 
lines  forming  the  "W"  should  be  straight.  However,  since 
the  upstream  waves  propagate  in  a  direction  inclined  by  © 

Si 

to  the  jet  axis,  these  dashed  lines  appear  with  a  small 
curvature. 

Periodic  oscillations  are  observed  in  the  correlation 
(4)  (t> 

plots  for  T  <T< T  in  Figs.  (4-15)  and  (4-16).  These 
oscillations  correspond  to  the  passage  period  2TT  /a»f. 
This  criterion  was  used  in  some  cases  to  define  peak  (4) 
in  the  correlation  functions. 
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The  prewhitened  correlation  functions  for  the  reso¬ 
nance  case  of  M  ■  .9  and  x^/d  "  4  are  shown  in  Fig. 
(4-17).  The  plots  show  a  similar  behavior  to  those  in 
Figs.  (4-1  A)  to  (4-16), 

Comparison  between  the  mathematical  model  and  the  op¬ 
timum  time  delays  from  the  measured  correlation  functions 
are  presented  in  Figs.  (4-18)  to  (4-20).  Eqns.  (4-17) 
and  (4-18)  can  be  written  in  the  nond imen sional  forms 


(3) 

q.g  Till, 


+  (m1C+  ™  e“)^r-  (4-i5) 


MKU  v  M  Kv 


<<0 

■°o  Titi  = 

A  X 


-  o>d  ^  +  cod.  ea)iii£, 

M  Kv  ax 


where  Ax  «  Xj  -  x^  is  the  separation  distance  along  the 

jet  axis  between  the  two  correlated  points  i  &  j.  The 

angle  ©  is  a  weak  function  of  the  plate  location  xn  [Fig. 
d  u 

_  O 

(4-12)],  therefore  an  average  value  of©  -30  is  used  to 
simplify  Eqns.  (4-25)  and  (4-26)  so  that  they  represent 
straight  lines.  Also,  as  shall  be  seen  in  the  following 
section  [Fig.  (4-21)],  the  convection  velocity  of  the 

downstream  wave  measured  from  the  optimum  tim*  delay 

u> 

7"  has  a  value  which  cor  r  e  ?po  r.'  ?•  to  Cj/U  *  ■  .62, 


I 


« 

i 

i 

i 
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while  the  upstream  wave  speed  is  Cg  *  a^. 

Fig.  (4-18)  shows  very  good  agreement  between  Eqn. 

<D 

(4-25)  and  the  measured  time  delay  7*  for  M  »  .8  and  .9. 

Even  at  low  subsonic  jet  speed  (M  ■  .5),  two  measured 

peaks  wf re  recognizable  in  the  correlation  functions  and 

are  in  good  agreement  with  the  model  in  Fig.  ( 4  —  1 8 b ) . 

These  points  indicate  the  existence  of  the  two  waves  even 

at  low  impinging  jet  speeds,  except  that  resonance  does 

not  occur.  Eqn.  (4-26)  is  also  in  excellent  agreement 

(4) 

with  measurements  of  7*  as  shown  in  Fig.  (4-19). 


The  magnitude  of  the  time  delay  of  the  extra  peak  in 
the  autocorrelation  function  from  Eqn.  (4-21)  can  be  ex¬ 
pressed  in  dimensionless  form  as 


(4.27) 


This  equation  is  plotted  in  Fig.  (4-20)  and  also  is  in 
good  agreement  with  the  measured  data. 


4.4  TWO  WAVES:  THE  DOWNSTREAM  AND  THE  UPSTREAM  TRAVELLING 

WAVES 

In  a  resonant  cas"  t“<:  nhase  velocity  at  the  reso¬ 
nance  frequency  car.  r-  I.- ted  from  the  cross-spectrum 
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between  two  point*  in  the  near  fields  For  point*  close  to 
the  jet  axis  (small  r) ,  end  point*  away  from  the  jet  axis 
(large  r)  •  Eqn .  (4-5),  after  some  manipulation  [see  Ap¬ 

pendix  A  and  Neuworth  (1973)],  takes  the  following  forms 
r/d  £  1 .0: 


-2b  K 


i  tfry*  -  M*-Xo>] 

e  +  (a-fe)  e  >  (4.28*) 


r/d  _>  1.0: 


/  L*rt  ([wt  +  K,.( 

«  jak  e  +  (b-o)  e  j  (4-jgb) 


where 


«=  ew.  j  (*-*)  J^r  (  rs  -  i-) 

The  first  term  on  the  right  hand  side  of  Eqn.  (4-2t) 
is  a  standing  wave.  The  second  tern  is  a  downstream  trav¬ 
elling  wave  in  Eqn.  (4-2|a)  and  an  upstream  travelling 
wave  in  Eqn.  (4-2|b).  Therefore,  at  resonance,  the  down¬ 
stream  phase  velocity  Cj(fr)  if,  ‘*t-  calculated  from  the 


t-.- '■  •  I  ’ — » — “ — rv 
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phase  difference  between  two  microphones  placed  near  the 
Jet  axis  (small  r) .  Also,  the  upstream  phase  velocity 
can  be  calculated  from  two  microphones  placed  in 
the  near  field  further  away  from  the  jet  axis  [see  Fig. 
(4-3}]  • 

The  phase  velocity  at  the  resonance  frequency  C(ff) 
and  the  broad-band  convection  velocity  C  for  the  two  waves 
are  plotted  in  Fig.  (4-21)  as  a  function  of  the  distance 
x  along  the  Jet  axis.  These  velocities  were  calculated  as 
follows : 


1.  The  downstream  phase  velocity  was  calcu¬ 
lated  from  Eqn.  (3-1)  using  the  phase  angle  ^(f^)  between 
two  points  closest  to  the  outer  edge  of  the  shear  layer 
(along  a  line  extending  from  the  norcle  lip  and  making  an 
angle  of  12*with  the  x-axis). 


2.  The  phase  velocity  of  the  upstream  wave  Cjffj.)  «t 
an  angle  of  8^  to  the  Jet  axis  was  calculated  from  the  re¬ 
lation 


C2  (fr)  * 


v  . .  ea 
% «« 


4  i  j  *  . y  )  1 1 
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The  angle  0^  is  given  by  Fig.  (4-12)  for  each  plate  loca¬ 
tion  x^/d. 

3.  The  downstream  and  upstream  convection  velocities 

Cj  and  were  calculated  from  the  correlations'  optimum 

(i)  JD 

time  delays  T  and  T  respectively. 

The  figure  indicates  that  the  upstream  wave  is  propa- 

A 

gating  in  a  direction  0  with  the  speed  of  sound. 
However,  the  downstream  wave  is  travelling  with  a  speed 
equal  to  .62  U.  This  is  the  convection  velocity  of  the 
large  scale  structures  as  measured  in  the  free  Jet  and  re¬ 
ported  in  Refs.  19  and  20.  Furthermore,  Neuworth  (1973) 
estimates  of  the  convecting  speed  of  large  eddies  in  a 
high  speed  impinging  jet,  obtained  from  a  frame  by  frame 
analysis  of  movie  film,  is  in  excellent  agreement  with  the 
present  date  as  shown  in  the  figure.  One  notices,  also, 
that  C ^ ( f r )  agrees  well  with  the  value  .62  U,  indicating 
that  large  coherent  eddies  play  the  main  role  in  the  reso¬ 
nance  phenomenon. 

The  conclusive  experimental  data  of  this  chapter  ex¬ 
hibits  two  waves  travelling  in  the  near  field  between  the 
nottle  and  the  plate.  Large  coherent  structures  are  con- 
vected  downstream  at  a  speed  equal  to  0.62  U.  Upstream 
waves  are  propagating  at  an  angle  0  to  the  jet  axis  with 


79 


the  speed  of  sound  Sg.  Furthermore,  the  simple  model  of 
Sec.  4.2  assumed  s  boundary  condition  of  zero  phase 
difference  between  the  two  waves  at  the  plate.  The  agree- 
■ent  between  the  measured  optimvim  time  delays  T  and 
T  and  the  model  implies  that  the  assumed  boundary  condi¬ 
tion  Is  indeed  correct.  We  conclude,  therefore,  that  the 
upstream  wave  Is  truly  a  reflection  from  the  plate  due  to 
the  Impinging  coherent  structures. 


CHAPTER  5 


THE  PHASE  LOCK  AND  THE  COLLECTIVE  INTERACTION 

The  establlshmen t  of  two  waves  achieved  In  Chapter  4 
is  not,  by  Itself,  sufficient  to  explain  the  resonance  ob¬ 
served  in  the  impinging  Jet.  In  order  to  have 
self-sustained  oscillations  and  to  close  the  feedback 
loop,  the  downstream  travelling  wave  and  the  upstream  pro¬ 
pagating  wave  should  be  phase  locked.  The  thin  shear 
layer  near  the  nozzle  exit  is  supposed  to  be  the  roost 
vulnerable  portion  of  the  jet  column  to  external  excita¬ 
tions.  Therefore,  the  shear  layer,  most  probably  near  the 
nozzle,  should  always  respond  to  the  upstream  forcing 
waves  thereby  closing  the  feedback  loop.  This  phase  lock 
has  never  been  established  experimentally  in  the  litera¬ 
ture:  it  is  the  main  feature  of  the  first  half  of  this 

chapter  . 

5.1  PHASF  LOCK  BETWEEN  THE  TWO  WAVES  AT  THE  NOZZLE  EXIT 

PLANE 

The  variations  of  the  rh?  >  -r.-'e  et  the  resonance 
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f  r  equency  ^>(  f  )  along  the  Jet  axis  for'both  waves  are  pre¬ 
sented  in  Fig.  (5-1).  The  figure  shows  the  case  M  »  .8, 

Xjj/d  ■  4  where  5506  Hz,  which  corresponds  to  (St)f  ■ 

.51.  The  variations  of^(fr)  along  the  downstream  travel-' 
ling  wave  were  determined  by  placing  a  series  of  micro¬ 
phones  close  to  the  shear  layer's  outer  edge  (r/d*1.25). 
The  phase  variations  from  the  plate  to  the  nozzle  exit 
plane  along  the  upstream  wave  were  calculated  from  micro¬ 
phones  at  r/d>2.5.  At  the  plate  the  two  waves  have  the 
same  phase  angle  (taken  equal  to  zero),  according  to  the 
boundary  condition  (4-9). 

The*7(fr)  variations  from  the  nozzle  lip  to  the  plate 
(at  rQ/d  "  1*0)  is  presented  in  the  figure  by  the  upper 

curve  and  the  lower  curve  repr esent s  ( f ^ )  variations  from 
the  plate  back  to  the  nozzle  exit  plane.  The  phase 
difference  between  the  two  waves  at  any  location  x  is, 
therefore,  equal  to  the  vertical  distance  between  the  two 
curves.  This  is  also  expressed  by  Eqn -  (4-11)  of  the 
model  in  Sec.  4.2.  The  first  term  represents  the  upper 
straight  line  with  a  slope  equal  to  kj  The  sec¬ 

ond  term  represents  the  lower  line  with  a  slope  propor¬ 
tional  to  1  /C2 •  The  most  important  result  of  the  figure 
is  that  the  phase  difference  between  the  two  waves  at  the 
nozzle  exit  (x  ■  0.)  is  a  multiple  of  27T.  This  indicates 
that  the  upstream  propagating  wave  is  in-phase  with  the 


downstream  travelling  wave  at  the  nosale  exit. 

The  conaiatency  of  the  above  aero  phaae  difference 
between  the  tvo  wave*  at  the  nossle  exit  vaa  examined  for 
different  plate  locations  and  is  ahovn  in  Fig.  (5-2). 
The  positions  of  the  four  near  field  microphones  used  are 
shown  In  Fig.  (4-3).  Two  microphones  were  placed  near 
the  Jet  (I  &  II)  to  weasur e  ^ ( f  )  variations  along  the 
downstream  travelling  wave.  The  other  two  microphones 
(HI4  XST  )  were  at  r/d  ■  2.7  to  measure  ^(fy)  variations 
over  the  upstream  propagating  wave.  For  each  plate  loca¬ 
tion  Xq,  the  phase  reference  point  on  the  plate,  which 
will  be  called  the  apparent  sound  source  "S",  is  deter- 
cined  as  follows  : 


1.  The  angle  between  upstream  propagation  and  the 
jet  axis  is  first  calculated  by  the  relation 


-I 

90  -  Co*. 


V  ;  (5-*> 


where  — - Sa5__l  i8  the  constant  phase  angle  slope  of  the 

d'f 

nondiapers ive  upstream  wave. 


2.  Using  ©  from  Eqn .  (5-1),  a  series  of  cross 

fl 

spectrum  functions  be t vec"  ZTT  and  different  surface  pres- 
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■ ur «  transducers  on  the  piste  were  evaluated.  The  appar¬ 
ent  sound  source  "S"  la  the  the  location  of  the  transducer 
which  best  aatlaflea  the  relation 

•in  - ±E>S_I — £»_  a  o0  .  (s-i) 

J/Jf 

other  words,  the  criteria  for  choosing  "S"  on  the  plate 
is  that  acoustic  propagation  from  S  to  m  In  a  direction 
making  to  the  jet  axis  be  given  by  Eqn .  (5-1), [and 

plotted  in  Fig.  (4-12)]. 

For  the  case  M  •  .9  in  Fig.  (5-2),  the  apparent 
sound  source  satisfying  Eqn.  (5-2)  was  found  to  be 

r#/d  -  1.0  for  3.25<  xQ/d  <7.5 

r#/d  -  1.5  for  xQ/d  <3.25 

The  percentage  deviation  from  satisfying  condition  (5-2) 
was  less  than  12X  for  transducers  within  one  nostle  diame¬ 
ter  of  the  above  values.  This  indicates  that  the  "source" 
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of  the  upstream  propagating  waves  emanates 


from  a  region 


on  the  plate,  rather  than  a  point 


The  phaae  variations  in  Fig.  (5-2)  for  all  Xg/d  show 
the  sane  characteristics  as  in  Fig.  (5-1).  The  phase 
difference  between  the  two  waves  at  the  nozzle  exit  rema¬ 
ins  as  an  Integer  muptiple  of  2Rfor  all  plate  locations. 
This  phase  lock  will  prove  to  be  Important  in  understand¬ 
ing  the.  resonance  frequency  stages,  the  topic  of  discus- 
sion  in  the  following  section. 

The  results  of  both  figures  indicate  that  the  2HN 
phase  difference  between  the  two  waves,  as  implied  by 
Powell's  (1961)  theory  on  edge  tones,  takes  place  at  the 
nozzle  lip.  This  closes  the  feedback  loop  which  allows 
for  a  self-sustained  system  of  oscillations.  Furthermore, 
the  results  of  the  conditional  sampling  technique  employed 
by  Petersen  (1978)  showed  that  the  pressure  measured  by 
microphones  at  the  outer  edge  of  the  shear  layer  is  90  out 
of  phase  with  the  cross  flow  velocity  component.  It  is 
well  known  that  this  phase  relation  is  also  the  same 
between  the  sound  wave  pressure  and  the  induced  particle 
velocity.  Therefore,  Bechert  and  Michel's  (1975)  in-phase 
assumption  for  the  cross  flow  velocity  of  a  shear  layer  at 
the  separation  point  of  a  semi-infinite  plate  and  the  vel¬ 
ocity  induced  by  external  acoustic  forcing  is  supported  by 
the  present  pressure  measurements. 
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5.2  RESONANCE  FREQUENCY  STAGES 


Although  the  phase  difference  between  the  two  waves 
at  the  nozzle  la  preaerved  for  different  plate  locatlona, 
the  reaonance  frequency  ff  changes.  Figs.  (5-3)  and 
(5-4)  show  the  variation  of  the  reaonance  Strouhal  number 
(St)^  •  f^d/U,  versus  the  plate  location  x^/d  for  M  ■  .8 
and  .9.  In  both  figures  (St)r  decreases  with  Increasing 
x^/d  until  It  reaches  a  minimum  value  [minimum  (St)  a»  .  5 
and  .33  for  M  *  .8  and  .9  respectively].  Beyond  the  mini¬ 
mum,  the  Strouhal  number  Jumps  to  a  higher  value,  then  de¬ 
creases  again  with  Increasing  xQ  and  the  cycle  repeats. 
The  frequency  Jumps  occur  over  a  small,  but  finite,  region 
of  plate  locations  such  as  4.8<Xg/d<5  and  6<Xp/d<6.5. 

At  H  «  .9,  a  strong  augmentation  of  energy  was  always 
observed  In  the  power  spectra  at  the  resonance  frequency 
[e.g.  Fig.  (3-22)].  Moreover,  this  energy  augmentation 
takes  place  at  a  single  resonance  frequency  most  of  the 
time,  except  when  the  plate  is  in  a  frequency  Jump  region 
as  Indicated  by  Fig.  (5-4).  However,  at  M  *  .8,  the  re¬ 
sonance  peaks  in  the  spectra  were  not  clearly  distinguish¬ 
able  for  x^/d  >3. 5  (also  screech  tones  were  barely  audi¬ 
ble).  This  can  be  explained  by  the  additional  frequency 
content  indicated  in  F ip.  (5-3';  clearly  the  energy  is 

86 


distributed  over  nore  then  one  frequency  band.  It  la  be¬ 
lieved  that  the  Impinging  jet  at  relatively  low  Mach 
numbers,  M<.8,  la  still  not  in  a  well  established  reso¬ 
nance  state.  This  Is  the  reason  for  choosing  the  case  of 
H  •  .9  to  study  most  of  the  details  of  the  resonance  phe¬ 
nomenon  in  the  present  investigation. 

The  resonance  frequency  stages  in  Figs.  (5-3)  and 
(5-4)  are  typical  of  self-sustained  oscillations  problems, 
such  as  flow  over  cavities  (Refs.  10,  13,  and  14)  and 
edge  tones  (Ref.  43).  Two  common  elements  are  used  in 
the  empirical  or  semi-eaplr leal  resonance  frequency  pred¬ 
iction  techniques.  A  phase  reference  is  assigned  and  an 
integer  is  used  to  fit  the  measured  frequency  at  different 
frequency  stages.  In  each  method,  the  question  of  when 
the  frequency  Jumps  occur  has  not  been  answered. 
Sarohia's  (1975)  prediction  of  experimentally  observed 
frequencies  is  based  on  a  prescribed  phase  difference 
equal  to  It  between  separation  and  impingement.  This  is 
only  one  branch  of  the  feedback  loop  which  means  that  his 
prediction  completely  neglects  upstream  acoustic  feedback. 
Rossiter's  (1964)  and  Block's  (1976)  models  include  the 
effect  of  an  acoustic  feedback.  However,  their  formulae 
require  an  empirical  constant,  which  is  in  effect  a  phase 
difference  between  the  feedback  acoustic  wave  and  the 
shear  layer  disturbances  near  the  upstream  separation,  as 
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vill  be  shown  later 


A  aero  phase  difference  between  the  two  waves  at  the 
noetic  exit,  as  established  by  the  present  data,  Implies 
that  the  number  of  wave. lengths  along  the  feedback  loop 
has  to  be  an  Integer*  The  feedback  loop  constitutes  a 
circuit  from  the  noztle  to  the  plate,  returning  again  back 
to  the  nozzle*  The  forward  branch  from  the  nozzle  to  the 
plate  is  due  to  large  scale  structures,  while  the  upstream 
acoustic  waves  form  the  rearward  branch  from  the  plate  to 
the  nozzle.  The  total  number  of  periods  at  any  time  over 
the  loop  Is  given  by 

N  =  —^2-  +  (5  -3 ) 

\«H  V.) 

where 

*.({.)=  ,  and 

tr 

V$r)A<*  6a  . 

*  tr 

are  the  wavelengths  of  the  two  waves  along  th  jet  axis. 
Eqn .  (5-3)  is,  in  f«.-*t,  Eqn  .  (4-10)  of  the  model,  with 
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■90  substituted  by  2KN 


Fig*  (5-4)  can  b*  raplottad  in  terms  of  N  veraua  the 
plate  location  xQ/d.  ThjL#  i§  ,hown  in  Fi8*  (5-5)  along 
with  the  conventional  plot  of  (St)f  versus  x^/d.  The  fig¬ 
ure  lndicatea  that  the  Integral  number  of  perioda  over  the 
feedback  loop  la  conatant  over  each  frequency  stage.  For 
Xg/dO,  near  field  measurements  were  not  poaaible  becauae 
of  interference  between  the  microphones,  ita  mountings, 
and  the  plate.  Eqn .  (5-3)  is  used  to  calculate  N  for 

xQ/d<3  using  the  values 

Cj(fr)  -  .62  II 

c2(fr)  -  «0 

a 

-  35  from  Fig.  (4-12) 

and  f ^  ia  measured  from  surface  pressure  transducers.  The 
calculated  values  of  N  follow  the  same  pattern  of  the 
measured  data  and  correspond  to  atage  N  ■  2  in  the  figure. 

The  frequency  stages  observed  in  flows  with 
ael f-austained  oscillations  has  constituted  a  pussl ing 
problem  for  some  time.  However,  with  the  knowledge  of  the 
convection  speeds  of  the  two  waves,  the  angle  of  propaga¬ 
tion  of  the  upstream  wave,  and  the  phase  lock, an  under¬ 
standing  of  the  frequency  stages  emerges  from  the  results 


of  Fig  *  (5-5) •  As  the  noxzl e- to-pl tte  separation  dis¬ 
tance  increasest  the  wavelength  of  both  waves  increases  to 
preserve  the  phase  lock  at  the  nozzle  exit.  This  results 
in  a  decrease  in  frequency  since  the  phase  velocity  is  al¬ 
most  constant.  A  minimum  frequency  that  corresponds  to  a 
long  wave  length  is  reached,  beyond  which  large  scale 
structures  cannot  maintain  their  coherence.  The  flow 
Jumps  to  a  higher  frequency  at  which  the  number  of  waves 
along  the  loop  Increases  by  one.  Then  the  frequency  de¬ 
creases  with  Increasing  separation  distance  and  the  cycle 
is  repeated. 


5.2.1  Prediction  of  Plate  Locations  at  the  Frequency  Jumps 


The  minimum  frequency  before  the  Jumps  corresponds  to 
(St>r  .33.  In  a  free  Jet,  large  scale  structures  are  ob¬ 
served  to  lose  their  coherence  beyond  the  end  of  the  po¬ 
tential  core.  The  passage  frequency  of  these  structures 
there  corresponds  to  the  above  minimum  Strouhal  number  be¬ 
fore  the  frequency  Jumps.  This  is  also  in  the  range  of 
the  Strouhal  number  of  the  preferred  mode  of  lsrge  scale 
structures  as  pointed  out  by  Crow  and  Champagne  (1971). 


The  minimum  frequency  can,  therefore,  be  used  with 
the  phase  lock  at  the  r  as  presented  by  Eqn .  (5-3) 


to  predict  plate  locat:  ■  'd  where  Jumps  occur.  Eqn. 


(5-3)  can  be  written  In  the  form 


x,  _  _J_  N 

J  (Sl)r  +  M  cm. 


(5-4q) 


oh*r€  m  Ij2j . 

Values  of  6  are  taken  from  Fig.  (4-12),  K  -  ,62  < nd 

A  V 

(St)  •  (St)  I  .  ■  .33  at  M  ■  .9.  Comparison  between  the 

plate  locations  at  which  frequency  jumps  occur,  as  pred¬ 
icted  by  Eqn .  (5-4a)  and  as  observed  from  the  measured 

data  of  Fig.  (5-4),  is  presented  in  Table  (5-1).  The 
predicted  values  in  the  second  column  is  in  excellent 
agreement  with  the  measured  values  in  the  third  column. 


— 

X 

3/d 

from  Eqn . (5-4a) 

from  Fig . (5-4 ) 

1 

1.25 

2 

2.  50 

2.4 

3 

3.  74 

3.75 

4 

4.99 

4.8  -  5.0 

5 

6.  24 

6.0  -  6.5 

6 

7.49 

7.5 

Table  (5-1)  Relative  distance  between  nozzle  exit 


and  plate  at  the  frequency  jumps. 


5.2.2  Jet  Behaviour  at  a  Frequency  Jump 

At  a  frequency  jump  more  than  one  reaonance  peak  ap¬ 
pears  In  the  power  spectrum  of  either  the  near  field  or 
surface  pressure  signals.  For  example,  the  spectra  of 
Fig. (4-3)  shows  two  peaks  representing  the  jump  at  x^/d-5 
in  Fig. (5-4). 

The  frequency  jump  phenomenon  is  studied  in 
Fig.  (5-6).  The  long  time  average  cross-spectrum  shows  two 
peaks  at  the  two  frequencies  that  correspond  to  the  jump 
at  xQ/d-2.35  .  When  the  spectrum  of  the  same  aignals  was 
calculated  over  short  time  intervals,  each  peak  intermit¬ 
tently  appeared  alone  .  This  indicates  that  during  a  fre¬ 
quency  jump  the  flow  switches  from  one  resonance  mode  to 
another,  and  this  happens  intermittently. 

The  mechanism  by  which  the  flow  field  switches  its 
frequency  during  a  frequency  Jump,  could  possibly  be  the 
same  as  the  one  observed  in  the  low  Reynolds  number  water 
jet.  There  the  frequency  at  which  vortex  rings  hit  the 
plate  intermittently  changed  as  the  flow  jumped  from  Pat¬ 
tern  1  (Fig.  (3-1)3  to  Pattern  2  (Fig,  (3-2)]  through 
"pairing"  . 
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5.2*3  The  Disappearance  of  Resonance 

The  resonance  disappeared  at  Xg/d  <  2  .  According  to 

Fig. (4-12),  the  upstream  propagating  waves  at  these  small 

plate  to  nozzle  distances  diverge  away  from  the  jet  axis. 

These  locations  correspond  to  a  minimum  frequency  stage  of 

N  .  -2  as  shown  in  Fig. (5-5). 
in  l  n 

At  large  plate  to  nozzle  distances  (x^  d>6), large 
scale  structures  would  lose  their  coherence  before  they 
Impinge  on  the  plate.  The  resonance.  consequently,  be¬ 
comes  weaker  and  then  disappears.  Furthermore,  an  oncom¬ 
ing  fully  developed  turbulent  flow  would  refract  a  wide 
band  acoustic  wave  away  from  the  jet  axis,  further  weaken¬ 
ing  the  feedback  to  the  nozzle  exit.  This  Is  shown  in  the 
plots  of  pressure  fluctuations  at  the  nozzle  exit  plane 
versus  Xg/d  in  Fig.  (5-7).  The  pressure  signals,  espe¬ 
cially  at  H  ,  are  dominated  by  the  upstream  acoustic  wave 
as  wa 8  shown  earlier  in  Ch .  4.  The  pressure  levels  are 

low  for  xQ/d>6.5  near  the  disappearance  of  resonance. 

5.2.4  Comparison  With  Cavity  Frequency 

Eqn.  (5-4a)  can  he  rewritten  in  the  form 
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l/ky  +  H  Co$>  9a 

N*  ijjj . 

The  above  equation  is  similar  to  Rosslter's  (1964) 
equation  for  cavity  frequency  which  takes  the  form 

lrx<>  _  n  -  r' 

u  »/Kv  +  M 

/ 

with  K  an  empirical  constant.  In  Eqn.  (5-4b)  the  empiri¬ 
cal  constant  is  absent  owing  to  the  zero  phase  difference 
between  the  two  waves  at  the  upstream  shear  layer  separa¬ 
tion  plane. 

In  Figs.  (5-8)  and  (5-9)  the  cavity  frequency  meas¬ 
ured  by  several  investigators  and  reported  by  Tam  and 
Block  (  1978)  is  compared  with  Eqn.  (5-4b)  .  Ir,  the  theo¬ 
retical  model  of  Tam  and  Block  (1978),  the  upstream  acous¬ 
tic  waves  reflect  from  both  the  bottom  and  the  upstream 
wall  of  the  cavity.  However,  they  argued  that  waves  re¬ 
flecting  off  the  upstream  wall  are  the  strongest  in  forc¬ 
ing  the  instability  of  the  shear  layer.  Therefore,  ©  is 

a 

taken  equal  to  zero  in  Eqn.  (5-4b)  since  the  present 


model  considers  the  feedback  to  occur  through  the  ambient 


region  (the  etetlonery  region  outeide  of  the  jet)  corres¬ 
ponding  to  the  Inside  of  the  cavity.  Eqn .  (5-4)  shows  a 
surprisingly  good  overall  agreement  with  the  measured 
data. 


If  the  shear  layer  of  the  Jet  is  approximated  by  a 
two-dimensional  one,  half  of  the  impinging  Jet  configura¬ 
tion  would  resemble  flow  over  a  cavity  of  infinite  depth. 
This  possibly  explains  why  Eqn.  ( 5  —  4 b )  appears  as  an 
upper  limit  to  the  measured  data  in  Figs.  (5-8)  and 
(5-9)  . 


5.3  THF,  COLLECTIVE  INTERACTION 


So  far  the  present  study  has  revealed  the  presence  of 
two  waves  in  the  near  field  of  the  impinging  jet.  The  up¬ 
stream  acoustic  waves  were  found  to  be  phase  locked  with 
the  downstream  travelling  hydrodynamic  waves  induced  by 
the  convected  coherent  eddies.  These  findings  close  the 
feedback  loop.  However,  a  few  questions  remain  to  be 
answered  concerning  the  coupling  between  the  upstream 
acoustic  waves  and  the  shear  layer  near  the  nozzle.  Do 
the  upstream  waves  force  the  initial  shear  layer  instabil¬ 
ity,  are  the  developed  instability  waves  different  from 
those  of  the  free  jet,  and  how  do  these  instability  waves 
evolve  downstream  into  large  coherent  structures  with  a 
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passage  frequency  equal  to  the  reaonance  frequency  ?  It 
waa  logical  then  to  aeaiure  the  lnatablllty  frequency  near 
the  nozzle  for  both  the  free  and  the  impinging  Jet. 

The  lnatablllty  frequency  of  the  free  jet  waa  meaa- 
ured  from  the  autocorrelation  function  of  a  pressure  sig¬ 
nal  ao  given  by  a  microphone  placed  at  the  nozzle  lip. 
This  frequency  is  plotted  for  different  jet  Mach  numbers 
in  Fig.  (5-12)  and  Is  discussed  later. 

For  the  impinging  Jet,  two  microphones  were  placed 
near  the  nozzle  exit.  The  raw  pressure  signals  are  shown 
in  Fig.  (5-10).  Fig.  (5-10a)  shows  signals  for  the 
non-resonant  case  of  M  *  .4  and  Xp/d  »  4.5.  A  high  fre¬ 
quency  signal  appears  immediately  downstream  of  the  nozzle 
exit  at  x/d  *  .13.  Farther  downstream  at  x/d  »  *92,  a  low 
frequency  signal  appears  with  the  high  frequency  signal 
superimposed.  The  low  frequency  signal  is  due  to  the 
evolving  large  scale  structures.  At  resonance  Fig. 
(5-l0b)  shows  that  for  M  *  . 9  and  x^/d  ■  4.5,  the  high 
frequency  components  are  superimposed  on  the  resonance 
frequency  signal  near  the  nozzle  exit.  At  the  downstream 
location  x/d  ■  1.31,  the  amplitudes  of  the  high  frequency 
signal  are  substantially  reduced. 

The  power  spectra  of  these  si;-  i ir.  shown  in  Fig. 
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(5-11)  and  ct'1  provide  inforaatlon  about  the  nature  of  the 
high  frequency  components  observed  in  the  pressure  sig¬ 
nals*  A  small  spectral  peak  at  high  frequency  appears  in 
the  spectrum  at  x/d  -  .13  in  Fig.  (5-lla).  The  same  peak 
becomes  clearer  in  the  spectrum  at  x/d  -  .92.  The  fre¬ 
quency  of  this  peak  is  exactly  the  instability  frequency 
measured  in  the  free  jet  (see  Fig.  (5-12)].  At  resonance 
Fig.  (5-llb)  shows  similar  high  frequency  spectral  peak, 
except  it  is  broader.  The  energy  contained  in  this  high 
frequency  peak  is  understandably  much  smaller  than  the  re¬ 
sonance  frequency  energy.  One  also  notices  that  unlike 
the  non-resonant  case  of  Fig.  (5-lla),  the  energy  at  high 
frequency  appears  to  vanish  at  the  downstream  location  x/d 
-  1.31. 

The  above  measurements  for  different  Jet  Mach  numbers 
are  plotted  in  terms  of  the  Strouhal  number  in  Fig. 
(5-12).  The  high  frequency  signals  for  both  the  free  jet 
and  the  impinging  jet  are  shown  to  follow  a  straight  line 
with  a  slope  equal  to  1/2.  This  is  in  agreement  with  the 
linear  instability  theory  of  Michalke  (1971).  The  data, 
thus,  indicate  that  the  initial  instability  frequency  of 
the  shear  layer  appears  unchanged  by  the  presence  of  the 
plate  or  the  feedback. 

The  resonance  frequency  is  also  shown  in  Fig. 
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(5-12).  The  vertical  bar*  uaed  in  the  plot  doe*  not  re¬ 
present  the  accuracy  of  tv^  measurements  but  rather  the 
range  of  resonance  frequencies  in  the  different  frequency 
stage,  as  it  appears  in  Figs.  (5-3)  and  (5-4).  The  fig¬ 
ure  shows  that  the  resonance  frequency  is  auch  lower  than 
the  instability  frequency  measured  near  the  nozzle  exit. 

The  ratio  (St).  /(St)  is  of  the  order  of,  or  even  larger 
m  r 

than  10.  The  large  difference  between  the  two  frequencies 
indicates  that  the  upstream  acoustic  waves  with  a  frequen¬ 
cy  f  cannot  be  phase  locked  with  the  hiRh  frequency  ini¬ 
tial  instability  waves  fin* 

The  mechanism  describing  the  drop  in  the  passage  fre¬ 
quency  of  large  eddies  in  a  two-dimensional  free  shear 
layer  was  first  reported  by  Winant  and  Browand  (1974). 
Pairing  between  vortex  structures  reduces  the  passage  fre¬ 
quency  by  a  factor  of  two  and  takes  place  over  about  one 
wavelength.  It  takes  at  least  3  to  4  pairings  then,  to 
reduce  the  frequency  by  a  factor  of  10.  Petersen  (1978) 
further  reported  that  the  passage  frequency  in  a  free  Jet 
decreases  linearly  with  downstream  distance  from  the  noz¬ 
zle.  This  is  not  true  in  the  impinging  jet  since  the  fre¬ 
quency  of  the  pressure  signals  of  Fig.  (5-10)  dropped 
sharply  over  a  short  distance  from  the  nozzle.  Therefore, 
conventional  pairing  cannot  be  adopted  to  explain  the 
sharp  drop  from  the  instability  frequency  to  the  resonance 
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one.  A  new  mechanism  is  needed  to  explain  the  situation 
near  the  nozzle. 

Ho  and  Huang  (1978)  observed  multiple  vortices  merg¬ 
ing  into  a  single  structure  when  a  low  frequency  driving 
force  is  applied  to  a  plane  shear  layer.  Fig.  (5-13) 
shows  one  of  their  dye  pictures.  The  flow  is  from  left  to 
right  with  the  higher  velocity  flow  at  the  bottom.  The 
low-frequency  forcing  displaces  the  vortices  from  their 
original  locations  and  leads  to  the  formation  of  an  array 
of  vortices  in  a  wavy  shear  layer.  Vortices  located  in  a 
portion  of  the  shear  layer  which  is  convex  as  seen  from 
the  low  speed  side  are  unstable  and  interact  together  to 
form  a  large  vortex. 

Ke  adopt  the  term  "collective  interaction"  to  des¬ 
cribe  the  phenomenon  of  multiple  merging  of  coherent 
structures.  The  characteristics  associated  with  this  phe¬ 
nomenon  are  a  sharp  drop  in  passage  frequency  and  a  rela¬ 
tively  large  shear  layer  growth.  Collective  interaction 
can  then  be  used  to  explain  the  frequency  drop  near  the 
nozzle  exit  in  the  impinging  jet. 

The  schematic  drawing  in  Fig.  (5-14)  can  assist  in 
explaining  the  collective  interaction  in  a  high  speed  imp¬ 
inging  jet.  Since  the  jet  is  r.N ;  ?  ymnetric  ,  only  the  upper 
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half  vi 11  be  diacuaaad.  The  ehear  layer  emerging  from  the 
nosele  le  puleating  due  to  periodic  forcing  from  the  up¬ 
stream  propagating  wave.  The  shear  layer  goes  through  cy¬ 
cles  of  divergence  (t  ■  0,  T  ,*•••)»  and  convergence  (t  * 

Tf /2 , 3T/2 ,♦••*) •  At  t  •  0  the  small  vortices  are  unst¬ 
able  with  respect  to  each  other,  and  tend  to  coalsce  Into 
a  large  vortical  structure  under  the  combined  effect  of 
their  induced  field  and  the  mean  shear.  At  t  »  Tr/2, 

these  merging  vortices  accelerate  in  the  curved  shear 
layer  and  are  subject  to  the  stabilising  effect  of 
stretching.  The  result,  as  shown  at  t  *  31^/2,  is  a  sharp 
drop  in  passage  frequency  from  the  instability  frequency 

f.  to  the  resonance  frequency  f  . 
in  *  r 

Looking  back  at  Fig.  < 5  —  1 0 b )  the  behavior  of  the 

pressure  signals  can  now  be  explained  in  the  following 
way:  the  high  frequency  signal  superimposed  on  the  low 

frequency  one  represent  the  instability  vortices  convected 
in  a  pulsating  shear  layer.  A  short  distance  downstream 
(x/d  •  1.31)  the  high  frequency  vortices  diminish  due  to 

the  collective  interaction.  Lau  et  al .  (1  972  )  have  re-  {* 

t 

lated  the  positive  peak  of  the  pressure  signal  to  the  pas¬ 
sage  of  the  valley  between  two  structures,  and  the  nega-  j 

! 

tive  peak  to  the  passage  of  the  bulge  of  the  structure.  1 

f  ' 

At  x/d  »  1.31  in  Fig.  (5— 10b) ,  the  amplitude  of  the  high 
frequency  waves  on  the  positive  per’s  of  the  low  frequency 
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waves  are  higher  than  those  on  the  negative  peaks.  These 
characteristics  further  confirm  the  collective  Interaction 
mechanism,  because  the  convex  part  of  the  shear  layer  to¬ 
wards  the  low  speed  side  (i.e.  the  bulge)  Is  less  stable. 
Therefore,  the  high  frequency  waves  on  the  bulge  decay 
faster  than  those  on  the  valley,  and  the  low  frequency  ne¬ 
gative  peeks  appear  smoother  than  the  positive  peaks  In 
the  figure. 

The  results  of  this  section  indicate  that  the  upstre¬ 
am  acoustic  wave  does  not  force  the  Initial  shear  layer 
instability  near  the  nozzle  exit.  The  forcing  acoustic 
waves  merely  drive  the  thin  shear  layer  In  phase.  Then, 
and  only  then,  does  the  wavy  shear  layer  act  on  the  high 
frequency  Initial  Instability  waves  and  collect  them  into 
large  coherent  structures  at  the  resonance  frequency. 

In  conclusion,  the  feedback  loop  between  the  nozzle 
and  the  plate  Is  presented'  in  the  flow  diagram  of  Fig. 
(5-15).  The  large  scale  coherent  structures  are  convected 
downstream  and  impinge  on  the  plate.  Waves  are  generated 
and  propagate  upstream  with  the  speed  of  sound.  These 
waves  drive  the  shear  layer  near  the  nozzle.  The  insta¬ 
bility  waves  undergo  collective  interaction,  while  the 
wavy  shear  layer  rolls  into  nev.*  larger  coherent  eddies  at 
the  resonance  frequency. 
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5.4  THE  RESPONSE  OF  THE  SHEAR  LAYER  TO  ACOUSTIC  EXCITATION 


Fl:OM  OUTSIDE  OR  INSIDE  THE  JET 

The  feedback  mechanism  described  in  the  previous  sec¬ 
tion  has  evolved  from  an  analysis  of  measurements  taken 
outside  the  Jet.  This  does  not  rule  out  the  possibility 
of  a  feedback  between  the  nozzle  and  the  plate  from  inside 
the  jet  column.  Wagner  (1971)  and  Neuworth  (1973)  adopted 
a  model  for  the  resonance  based  on  waves  travelling  in  the 
jet  core,  supported  by  visualizations  of  standing  waves  in 
the  jet  column.  The  present  investigation  does  not  rule 
out  a  feedback  within  the  jet.  In  fact,  a  phase  lock  of 
Ttbetween  the  acoustic  waves  propagating  inside  the  jet  and 
outside  the  jet  would  be  more  powerful  in  driving  the 
shear  layer  near  the  nozzle.  Unfortunately,  measurements 
inside  a  high  subsonic  speed  jet  are  difficult  and  unreli¬ 
able.  Hussain  and  Zaman  (1977)  have  shown  that  the  probe 
itself  can  trigger  an  upstream  instability  mode  when  in¬ 
serted  in  a  free  shear  layer.  Neuworth  (1973)  argued  that 
feedback  could  even  occur  due  to  a  small  perturbing  body 
as  long  as  it  has  a  stagnation  surface  within  the  flow. 

In  this  section  a  theoretical  approach  is  used  to  de¬ 
termine  the  relative  importance  of  acoustic  feedback  in- 
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•id*  th«  jat.  Tht  *ff«ctlvn««a  of  acoustic  waves  in 
triggering  the  shear  layer  instability  near  the  noerle 
from  outside  the  jet  (casa-A)  is  compared  to  that  from  in¬ 
side  the  jet  core  (casa-B) .  The  two  cases  (A  and  B)  are 
shown  schematically  in  Fig.  (5-16). 

For  simplicity,  consider  plane  acoustic  waves  propa¬ 
gating  in  a  two- dimenaional  flow  field  on  both  sides  of  n 
zero  thickness  shear  layer.  A  aero  thickness  shear  layer 
is  Justified  since  the  acoustic  wavelength  is  much  larger 
than  the  anticipated  shear  layer  thickness.  Quantitative 
correctic..  factors  due  to  finite  shear  layer  thickness  and 
cylindrical  flow  field  are  reported  by  Amiet  (1977). 
These  corrections,  however,  do  not  alter  the  physical  pic¬ 
ture  of  the  present  analysis. 

The  approach  is  based  on  the  theoretical  work  of 
Ribner  (1957).  The  Induced  cross  flow  velocity  at  the  in¬ 
terface  V  due  to  the  incident  acoustic  waves  pressure  p^ , 
is  a  measure  of  the  forcing  of  shear  layer  instability. 
For  the  same  incident  acoustic  wave  pressure,  the  ratio  of 


cross  flow  velocities  in  the  two  cases  is 


where 
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It!  And  0t  are  the  transmitted  Acoustic  wive  Ampli¬ 
tude  And  Angle  to  the  interfsce.  The  Amplitude  It!  tskes 
Lhe  form 
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Umo  Sin.  2&i  -+»  sin.  £ 


The  rAtio  V  /V  is  plotted  for  a  range  of  incidence 

•  a/b 

Angles  0  <  ©i  <90  in  Fig.  (5-16).  This  is  the  rsnge  of 

concern  here  since  it  represents  only  upstream  Acoustic 

propsgsting  waves.  At  normal  incidence  (0.  »  90  )  V  /V 

l  CA  lB 

is  unity.  As  0^^  decreases,  the  rAtio  increases  sharply  to 
infinity  when  ©  1  is  equal  to  a  critical  angle  ©c.  Below 
,  a  complete  sound  reflection  occurs  in  case-B  and  the 
shear  layer  acts  as  a  sound  hArd  boundary.  For  a  typical 
resonant  case  of  It  •  .9,  the  critical  angle  is  as  high  as 
63  .  Thus,  for  a  considerable  range  of  incidence  angle, 
the  shear  layer  is  far  more  vulnerable  to  upstream  acous¬ 
tic  waves  from  outside  the  jet,  especially  at  high  Mach 
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The  above  simple  results  suggest  that  the  present 
measurements  outside  the  jet  are  Indeed  sufficient  to  pro¬ 
vide  a  complete  picture  of  the  feedback  mechanism. 


CHAPTER  6 


FAR  FIELD  NOISE  OF  THE  IMPINGING  JET 


The  far  field  noise  of  the  in pinging  Jet  for  both  re¬ 
sonant  and  nonresonant  cases  is  investigated  in  this 
chapter.  Powell  (1960)  suggested  that  the  rigid  surface 
plays  a  passive  role  in  the  noise  radiation.  The  noise  is 
generated  by  the  flow  and  is  being  reflected  by  the  plate. 
A  mirror  image  of  the  flow  could  replace  the  plate  in  gen¬ 
erating  far  field  noise.  Therefore,  in  his  model,  the 
flow  and  the  plate  have  equal  shares  in  the  noise  radia¬ 
tion.  Preisser  (1979)  measured  the  noise  spectra  of  the 
Impinging  jet.  The  spectral  peaks  at  (St)  *  .3,  were 
higher  in  magnitude  than  those  in  the  spectra  of  a  free 
jet.  He  suggested  that  the  presence  of  the  plate  had  en¬ 
hanced  the  noise  generated  by  the  large  scale  structures. 
Most  of  the  noise  is  reported  to  be  generated  in  the  stag¬ 
nation  region  near  the  plate.  However,  Preisser  could  not 
determine  whether  the  noise  was  generated  in  the  jet  and 

reflected  by  the  plate  or  whether  It  was  generated  due  to 
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the  impinging  of  the  Jet  on  the  plate.  Pan  (1975)  pro- 
poaad  a  theoretical  method  baaed  on  correlation  techniquea 

1 

I 

to  aaparat e  the  turbulent  flow  contributlona  to  the  far  • 
fiald  nolae  from  those  due  to  surface  pressure  fluctua¬ 
tions.  He  expressed  the  need  for  near  field  measurements 
as  well  as  surface  pressure  measurements  to  determine  the 
exact  role  of  the  plate  in  the  generation  of  noise. 

The  present  far  field  measurements  are  coupled  with 
near  field  and  surface  pressure  measurements  to  determine 
the  mechanism  of  noise  generation  and  its  radiation  path. 

The  far  field  measurements  are  taken  from  a  microphone  at 

e 

©■  89  i.e.,  approximately  in  a  plane  perpendicular  to  the 
Jet  axis  at  the  nozzle  exit.  Shielding  of  the  plate  pre¬ 
vented  reliable  measurements  at  smaller  angles,  while  that 
of  the  nozzle  assembly  prevented  large  angle  measurements. 

Fig. (6-1)  shows  the  coordinate  system  that  will  be 
used.  X  is  used  here  as  a  composite  coordinate  that  ex¬ 
tends  from  the  nozzle  along  the  Jet  axis  to  the  stagnation 
point  on  the  plate,  and  then  radially  on  the  plate.  This 
is  approximately  the  path  of  Induced  disturbances  in  the 
near  field  and  on  the  plate  of  an  impinging  vortical  ring 
[Ref.  33].  Hence ,  1  sa t i f ie s  the  following  relations: 

^  -  x  in  the  near  field> 

\  -  Xq  +  ( rQ  -  r  1  on  the  plate. 


In  the  figure  C  is  the  distance  between  any  point  in  the 
near  field,  or  on  the  plate,  and  the  far  field  microphone. 

6.1  THE  ROLE  OF  LARGE  SCALE  STRUCTURES  IN  NOISE  GENERATION 

2 

The  overall  far  field  sound  pressure  level  p'  ,  meas¬ 
ured  in  decibels,  are  shown  in  Fig.  (6-2)  as  a  function 

/  2 

of  the  plate  location  Xp/d  for  M  ■  .9.  The  value  of  p 
fluctuates  around  135dB  and  then  falls  off  to  llOdB  for 
Xg/d>6.5.  The  variations  of  pf  ^  with  Jet  Mach  number  M 
are  shown  in  Fig.  (6-3)  for  x^/d  -  4.5.  p/  increases 

g 

proportionally  to  M  as  in  a  free  Jet.  The  sound  pressure 
level  for  a  free  jet  at  M-0.8  is  included  and  is  seen  to 
approximately  equal  that  of  the  impinging  jet.  This  is  in 
agreement  with  measurements  by  Preisser  (1979),  indicating 
that  at  directions  perpendicular  to  the  jet  axis  the  noise 
levels  are  almost  the  same  for  both  the  free  and  the  imp¬ 
inging  jets.  However,  at  large  angles,  Preisser's  meas¬ 
urements  do  show  higher  noise  levels  for  the  case  of  an 
impinging  jet  indicating  that  the  presence  of  the  plate  is 
felt  more  in  the  far  field  at  large  angles. 

The  mechanism  of  noise  radiation  is  investigated  here 
using  the  correlation  functions  between  far  field,  near 
field,  and  surface  pressure  measurements.  Fig.  (6-4) 


•hows  the  correlations  between  the  far  field  microphone 

and  three  near  field  microphones.  Also,  the  correlation 

between  the  far  field  microphone  and  a  pressure  transducer 

on  the  plate  at  r^/d  ■  1.0  la  presented.  The  time  delay 

In  the  figure  starts  from  sn  offset  delay  T.  ■  5.625  msec. 

a 

In  order  to  present  only  the  doainsnt  features  of  the 

correlation  functions.  A  negstlve  optimum  peak  appears  In 

the  correlations  with  the  near  field  microphones,  while  a 

positive  one  appears  In  the  correlation  with  the  transduc- 
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er  on  the  plate.  The  time  delay  T  of  the  peaks  of  the 
correlations  with  the  nesr  field  satisfies  the  lnea.uallty 

(i)  X 

Tj,rr  >  *j/o ©  «..*) 

The  above  inequality  indicates  that  the  time  delay  between 
any  point  in  th  near  field  and  the  far  field  is  longer 
than  the  time  it  takes  for  direct  acoustic  radiation 
between  the  two.  Therefore,  the  possibility  of  an  acous¬ 
tic  path  from  the  near  field  directly  to  the  far  field  is 
exclud ed  . 

For  different  jet  Mach  numbers  and  plate  locations 
Xg/dv  the  following  relations  were  always  true: 
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The  terms  between  brackets  are  the  optimum  time  dela¬ 
ys  of  peak  (1)  in  the  correlations  of  Fig.  (4-16)  corres¬ 
ponding  to  the  downstream  convection  of  large  scale  struc¬ 
tures.  Eqn.  (6-3),  therefore,  indicates  that  the  time 
delay  between  any  point  in  the  near  field  and  the  far 
field  is  equal  to  the  summation  of  near  field  time  delays 
up  to  point  2  on  the  plate,  plus  the  time  delay  between 
point  2  and  the  far  field.  Moreover  Eqn.  (6-4)  shows 
that  the  time  delay  between  point  2  (in  the  region  where 
large  scale  structures  impinge  on  the  plate)  and  the  far 
field  corresponds  to  direct  acoustic  radiation. 

According  to  Eqns.  (6-3)  the  role  of  the  plate  as  a 
reflector  of  sound  generated  by  the  flow  as  indicated  by 
Powell's  (1960)  model  is  incorrect.  The  time  delay  terms 
in  brackets  in  Eqn.  (6-3)  should  correspond  to  a  convec¬ 
tion  velocity  that  is  equal  to  the  speed  of  sound  in  order 
for  Powell's  model  to  be  correct.  However,  these  time  de¬ 
lays  correspond  to  .62  U,  far  below  the  speed  of  sound  a^ 
(e.g.,  at  M  *  .8,  .62U*.5aQ). 

The  results  of  the  experimental  data  as  indicated  in 
Eqns.  (6-3)  and  (6-4)  explain  the  noise  generation  and 
its  radiation  path.  The  large  scale  structures  interact¬ 
ing  with  the  plate  pi  3v  -  i-.-.rmnt  role  in  the  genera- 
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tlon  of  noise.  The  nolae  does  not  radiate  directly  from 
the  downstream  convected  large  coherent  structures,  but 
rather  as  these  structures  impinge  on  the  plate.  That  is 
why  the  surface  and  far  field  measurements  of  Preisser 
(1979)  indicated  a  strong  apparent  noise-producing  region 
within  one  jet  diameter  of  the  stagnation  point  on  the 
plate . 


Petersen  (unpublished)  cross  correlated  signals  from 
microphones  in  the  near  and  far  fields  of  a  free  Jet.  The 
time  delays  of  the  correlation  peaks  Indicated  that  pres¬ 
sure  fluctuations  propagate  down  the  Jet  column  and  radi¬ 
ate  sound  from  a  compact  acoustic  source  region  near  the 
end  of  the  potential  core.  However,  the  propagation  speed 
of  the  pressure  fluctuations  in  the  near  field  was  high 
(0.88U).  Since  Michalke's  (1971)  theory  indicates  that 
for  small  amplitude  disturbances  long  waves  travel  faster 
than  short  waves,  Petersen  suggested  the  possibility  that 
hydrodynamic  disturbances  relevant  to  the  noise  production 
may  be  fast  moving  long  waves.  Except  for  this  question¬ 
able  high  convection  velocity,  it  is  interesting  to  note 
that  his  model  for  the  noise  radiated  from  a  free  jet  is 
in  a  sense  similar  to  ours.  The  noise  seems  to  radiate 
from  near  the  end  of  the  potential  core  where  pairing 
between  large  structures  usually  takes  place  and  somehow 
resembles,  the'  impingement  of  large  structures  on  the 
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plate  In  our  case. 

The  peak  (1)  in  the  correlation  functions  between  the 
near  and  far  field  microphones  la  shown  In  the  figure  to 
be  negative,  while  it  is  positive  In  the  correlation 
between  the  far  field  and  the  surface  signals.  With 
reference  to  Fig.  (6-4),  the  following  explanation  of 
this  phenomenon  Is  offered.  Consider  a  frame  of  reference 
moving  downstream  with  the  large  scale  structures.  A 
large  eddy  at  I,  for  example,  would  be  stationary  with  a 
flow  convected  over  it  with  velocity  (from  right  to 

left).  The  flow  accelerates  as  It  moves  over  the  top  of 
the  eddy,  while  its  pressure  drops  according  to  Bernoulli 
equation.  A  microphone  at  I  would  then  record  a  negative 
pressure  disturbance.  This  is  supported  by  the  near  field 
measurements  of  Lau.et  al  (1972),  who  showed  that  the  ne¬ 
gative  pressure  disturbance  corresponds  to  the  passage  of 
a  "bulge"  in  the  shear  layer.  When  this  eddy  Impinges 
downstream  on  the  plate,  a  positive  pressure  disturbance 
1 8  recorded  at  2  and  radiated  to  the  far  field.  The  mi¬ 
crophone  in  the  far  field  measures  the  positive  distur¬ 
bance  and,  accordingly,  a  positive  peak  appears  in  the 
correlation  with  the  plate  Rj  Fp»  On  the  other  hand,  the 
same  positive  disturbance  at  FF  has  been  measured  earlier 
as  a  negative  one  at  I,  resulting  in  a  negative  peak  in 

RI,FF* 
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A  second  peak  appears  In  the  figure  and  la  denoted  by 
(11)*  In  the  correlations  between  FF  and  other  more  dis¬ 
tant  downstream  locations  (II  and  III),  peaks  (1)  and  (11) 
move  towards  each  other  and  ultimately  come  Into  coince- 
denceclde  In  the  correlation  with  the  plate  Rj  pp*  The 
time  delay  of  the  second  peak  f  was  found  to  satisfy 
the  following  relation 


r<0 
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The  right  hand  side  of  Eqn  .  (6-5)  is  the  time  delay  of 
peak  (2)  in  Fig.  (4-16),  which  corresponds  to  the  upstre¬ 
am  acoustic  propagation.  Peak  (ii)  in  the  correlation 
functions  of  Fig.  (6-4)  can,  therefore,  be  explained  as 
follows:  a s  the  large  scale  structures  impinge  on  the 
plate,  acoustic  waves  are  not  only  radiated  ti.'  the  far 
field  FF,  but  also  to  I  on  account  of  the  near  fir'd  up¬ 
stream  waves.  While  a  peak  (i)  in  RT  __  i*»  due  to  the 
path  from  I  to  2  to  FF,  peak  (ii)  corresponds  to  the  same 
path,  minus  the  upstream  propagation  from  2  to  I.  So  as  I 
moves  downstream  (to  II  and  III),  the  difference  between 
the  time  delays  of  the  two  peaks  decreases  due  to  the 
shorter  upstream  propagation  distance.  When  I  reaches  the 
plate,  the  two  peaks  merge  into  one.  It  should  be  menti¬ 
oned  here  that  peak  (ii)  does  not  represent  direct  radia- 


tlon  from  the  large  scale  structures  to  the  far  field. 

The  radiation  mechanism  corresponding  to  the  measured 
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r  and  T  will  be  further  confirmed  later  by  comparison 
with  a  model  based  on  Eqns.  (6-3),  (6-4)  and  (6-5). 


6.1.1  Major  Noise-Production  Region  on  the  Plate 

The  correlation  functions  between  the  far  field  mi¬ 
crophone  and  the  pressure  transducers  on  the  plate  are 
plotted  in  Fig.  (6-5).  Correlations  with  points  in  the 
developed  wall  jet  region  are  too  weak  and,  therefore,  are 
not  presented.  The  plots  are  presented  for  time  delays 
larger  than  5.781  msec.,  since  no  significant  correlation 
levels  are  observed  before  that.  The  level  of  the  corre¬ 
lation  peak  decreases,  and  the  width  slightly  increases 
with  distance  from  the  stagnation  point  on  the  plate.  The 
increase  in  width  is  indicative  of  the  stretching  of  large 
vortices  during  and  after  impingement.  Based  on  the  dis¬ 
tance  to  the  far  field  microphone  ,  and  the  speed  of 
sound  8q,  the  time  delays  of  these  peaks  do  not  correspond 
to  direct  acoustic  radiation  except  for  point  2. 
Furthermore,  the  optimum  time  delay  at  any  point  on  the 
plate  satisfies  the  relation 


T 


C 


(on  hht  plait)  (6-  it) 


where  C  ie  the  eddy  convection  velocity  elong  the  plate, 
presented  earlier  In  Fig.  (3-26).  As  seen  from  Eqn. 
(6-6),  the  nolae  radiates  from  the  impinging  region  around 
r^/d  ■  .  on  the  plate  (represented  by  point  2).  At  points 
outside  that  region,  the  time  delay  differs  by  an  amount 
corresponding  to  the  convection  of  large  eddies  from  point 
2  to  that  point. 

The  large  correlation  coefficient  in  Rq  ^(T)  of  Fig. 
(6-5)  is  not  due  only  to  high  surface  pressure  fluctua¬ 
tions  at  Tp/d  •  0,  since  the  same  was  observed  for  small 
Xp/d,  where  surface  pressure  fluctuations  at  the  stagna¬ 
tion  point  were  low.  In  order  to  explain  the  high  level 

peak  in  R.  __(T),  one  should  take  into  consideration  that 
o ,  F  F 

large  cross-correlation  coefficients  not  only  imply 
phase-related  signals,  but  also  overlapping  spectral  con¬ 
tents.  These  results  further  support  our  earlier  observa¬ 
tion  [see  Fig.  (3-21)]  that  the  induced  signal  at  rp*0 
depicts  the  gross  features  of  the  impinging  vortical 
structures.  Therefore,  it  is  not  surprising  that  this 
particular  point  is  strongly  correlated  with  the  far 
field.  We  conclude  from  Fig.  (6.5)  and  Eqn.  (6-6)  that 
the  inner  impinging  region  is  the  major  noise  source  on 
the  plate,  in  agreement  with  Preisser  (1979). 
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6.1*2  A  Model  for  the  Noise  Generated  by  Large  Coherent 


Impinging  Structures 

The  behavior  of  the  correlations  of  the  far  field  mi¬ 
crophone  with  both  the  near  field  microphones  and  surface 
pressure  transducers  can  be  explained  by  the  sketch  of 
Fig.  (6-6).  The  noise  is  mainly  generated  by  the  im¬ 
pingement  of  the  large  scale  structures.  The  impinging 
region  on  the  plate  is  represented  by  the  apparent  sound 
source  "S".  "S"  has  been  identified  for  different  plate 

locations  in  Chapter  5,  and  was  found  in  most  cases  to  lie 
near  rQ/d  «  1.  After  impingement,  the  large  eddies  con- 
vect  along  the  plate,  while  a  wave  propagates  upstream  in 
the  near  field. 

Eqns.  (6-3),  (6-4)  and  (6-6)  can  be  reduced  to  the 
nond imensional  equation 

IA.,1  -I  -  1—  1  ,  (*-T) 

where  the  composite  coordinate  of  Eqn  .  (6-1)  is  used. 

Implicit  in  Eqn.  (6-7)  Is  the  assumption  of  Equal  convec¬ 
tion  speeds  in  the  near  fit  Id  and  along  the  plate.  This 
is  Justified  since  the  i  rastion  is  to  be  examined  using 


surface  measurements  near  the  Impinging  region  on  the 
plate  where  the  eddy  convection  velocity  has  the  average 
value  0  *  62U  (see  Fig.  (3-26)]. 


The  unit  coefficient  appearing  aa  the  first  term  in 
the  right  hand  side  of  Eqn.  (6-7)  repreaents  the  direct 
acoustic  radiation  from  "S".  The  second  term  is  due  to 
the  convection  of  large  eddies  to  or  from  "S". 

W 

The  time  delay  of  the  second  peak  T  (Fig*  (6-4)] 
can  be  derived  from  Eqn.  (6-5)  and  put  in  the  nondimen- 
sional  form 


(it) 

T.  °9  s  I  -  — 1-, * _  CoS.  ©a  (near  field  only).(6-g) 
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The  second  term  on  the  right  hand  side  is  due  to  the  near 
field  upstream  wave  propagating  at  an  angle  to  the  Jet 
axis  . 


Eqn s .  (6-7)  and  (6-8)  are  compared  with  the  measured 
data  in  Fig.  (6-7).  The  cases  of  x ^ / d  •  5.5  and  7  at  H  « 
.8  are  shown.  A  value  of  ©^  *  27*  is  used  for  both  cases 
as  given  by  Fig.  (4-12).  T  c  apparent  sound  source  is  at 
r^/d  «  1.0  on  the  plat  .  •  fured  data  points  at 


015,  0, 


(  h  **  •*  )  “  .015,  0,  -.  0055  and  -.012  correspond  to  corre- 

latlon  functions  with  points  on  the  piste.  It  is  evident 
that  the  measured  data  are  in  excellent  agreement  vlth 
Eqns.  (6-7)  and  (6-8). 


The  results  discussed  so  far  indicate  that  the  noise 
is  generated  by  the  impingement  of  large  scale  structures 
on  the  plate.  Accordingly,  the  power  spectrum  of  the  far 
field  signal  is  expected  to  peak  at  a  low  frequency  band 
that  corresponds  to  the  passage  frequency  of  these  struc¬ 
tures.  The  far  field  signals  were  digitized  at  a  very 
high  sampling  rate  to  retain  any  possible  high  frequency 
content.  A  presentation  of  a  "surprising  behavior"  of  the 
power  spectra  is  given  in  the  next  section. 


6,2  HIGH-FREQUENCY  FAR  FIELD  NOISE  RADIATION 

Fig. (6-8)  shows  two  normalized  power  spectra  G(f)  for 
Xp/d  *  4.5.  The  spectra  were  evaluated  from  data  that 
were  digitized  at  a  fast  sampling  rate  of  4.7  microsec. 
The  non-resonant  case  at  H-.4  is  plotted  in  Fig.(6-8a). 
The  spectrum  has  a  peak  at  f-2000  Hz  corresponding  to 
St«*.33.  This  value  falls  within  the  frequency  range  of 
the  large  scale  structures,  a  fact  which  supports  the 
noise  radiation  mechanism  discussed  earlier.  However,  a 
considerable  amount  of  ercr?  ■  •  «■  or.tr.ined  at  high  fre- 
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quencies,  as  indicated  by  the  second  peak,  around  25  RHz. 

A  sharp  spectral  peak  also  appears  at  f*17.7  RHz  this  . 
being  exactly  the  instability  frequency  measured  at  the 
nozzle  lip  (Sec*  5.3). 

A  similar  behaviour  is  shown  in  Fig.(6-8b)  for  M*.9. 
The  high  frequency  peak  is  now  much  smaller  than  the  low 
frequency  peak.  This  is  not  unexpected  since  in  this  re¬ 
sonant  case,  most  of  the  energy  is  concentrated  at  the  re¬ 
sonance  frequency  t  (St)  *.35  ),  a  point  which  will  be 
discussed  later.  Note  that  a  spectral  peak  also  appears 
at  the  resonance  frequency's  first  harmonic  ( 2  f  ^ )  • 

It  wa 8  shown  earlier  that  the  frequency  of  the  waves 
in  the  shear  layer  drops  sharply  over  a  short  distance 
from  the  nozzle  due  to  the  collective  interaction.  As  a 
result,  large  structures  convect  downstream  and  impinge  on 
the  plate,  radiating  sound  to  the  far  field.  The  follow¬ 
ing  questions  may  then  be  raised:  From  which  location 
does  the  high  frequency  noise  in  the  measured  spectra  ori¬ 
ginate  ?  How  does  it  radiate  ?  It  seems  reasonable  to 
look  for  high  frequency  sound  radiators  close  to  the  noz¬ 
zle  exit,  and  to  expect  a  different  sound  path  from  the 
one  due  to  the  low-frequency  large-scale  structures  dis¬ 
cussed  in  the  previous  section.  Correlations  between  the 
far  field  microphone  and  near  field  microphones  placed 
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along  the  outer  edge  of  the  shear  layer  and  near  the  noz¬ 
zle  exit  were  used  to  investigate  the  high  frequency  noise 
radiation. 

An  example  of  the  correlation  with  a  near  field  mi¬ 
crophone  at  x/d*.92  for  the  case  of  M«.4  and  Xg/d-A.S  is 
plotted  in  part  (a)  of  Fig(6-9)  *  A  broad  negative  peak 
similar  to  peak  (i)  in  the  correlation  functions  of 
Fig. (6-4)  is  shown.  The  time  delay  of  this  peak  satisfies 
Eqn.(6-7),  thereby  indicating  that  it  is  associated  with 
the  noise  generated  by  the  impingement  of  large  scale 
structures  on  the  plate.  Peak  (ii)  does  not  appear, 
however,  because  of  the  weak  feedback  from  the  plate  for 
such  a  low  jet  speed.  A  peak  denoted  by  (ill)  emerges  at 
an  earlier  time  delay.  The  narrow  width  of  this  peak  in¬ 
dicates  a  correlation  between  high  frequency  events. 

To  ensure  that  the  above  correlation  peak  corresponds 
to  the  radiated  high  frequency  noise,  the  pressure  signals 
were  passed  through  a  digital  high-pass  filter.  The  low 
frequency  cutoff  frequency  of  the  filter  was  chosen  to  be 
the  8  KHz  frequency  of  the  valley  between  the  two  spectral 
peaks  in  Fig.  (6-8).  The  correlation  between  the  filtered 
signals  is  shown  in  part  (b)  of  Fig. (6-9).  Although  the 
broad  peak(i)  disappeared,  peak(iii)  due  to  the  high  fre¬ 
quency  components  of  the  sipnals  is  retained. 
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Peak  (iii)  could  not  be  identified  in  correlations 
with  farther  downstream  microphones  in  the  near  field  . 
This  Indicates  that  high  frequency  noise  is  generated  ma¬ 
inly  near  the  nozzle  exit,  before  collective  interaction 
takes  place*  The  region  of  high  frequency  fluctuations  in 
the  shear  layer  gets  smaller  with  strong  feedback  at  high 
jet  speeds  t  Fig. (5-11)  ].  Hence  a  relatively  small  high 
frequency  spectral  peak  is  observed  at  M«.9  in  comparison 
to  the  case  M*.  4  of  Fig.  (6-8). 

The  noise  radiation  path  can  be  determined  from  the 
time  delay  of  peak  (iii).  This  is  plotted  in  Fig.  (6-10) 

W(H»  W 

where  t  is  the  distance 
between  the  near  field  and  the  far  field  microphones. 
Cases  M  *  .3,  .4,  and  .5  are  shown.  The  data  collapse 
onto  a  straight  line  at  *  I -0 ,  indicating  direct  high 
frequency  acoustic  radiation  to  the  far  field. 

Note  that  at  H  ■  . 5,  only  radiation  from  I  (next  to 
the  nozzle  exit)  is  detected  and  shown  in  Fig.  (6-8)  This 
corresponds  to  the  fact  that  high  frequency  radiation  is 
restricted  to  shorter  distances  from  the  nozzle  exit  when 
the  jet  speed  is  increased.  At  higher  subsonic  jet  speeds 
resonance  dominated  the  correlations.  When  the  prewhiten¬ 
ing  technique  was  used,  a  r*' sk  a1  ways  appeared  at  a  time 


delay  correpondlng  to  direct  noise  radiation*  However, 
the  level  of  this  peak  was  low  and  within  the  correlation 
levels  of  the  noise  in  the  signals.  As  expected,  this  is 
due  to  the  low  energy  content  of  the  high  frequency  spec¬ 
tral  peak  as  shown  in  Fig.  (6-8b)  for  M  ■  .9. 
Consequently,  data  points  at  high  Mach  numbers  were  not 
plotted  in  Fig.  (6-10). 

A  literature  survay  did  not  reveal  any  previous  evi¬ 
dence  for  the  presence  of  two  peaks  in  the  far  field  spec¬ 
tra  for  either  a  free  jet  or  an  impinging  jet.  Preisser's 
(19/9)  impinging  jet  has  a  thick  initial  shear  layer  due 
to  the  use  of  a  long  nozzle.  Since  the  frequency  of  the 
instability  waves  scales  with  the  shear  layer  thickness  at 
the  nozzle  exit,  high  frequency  noise  radiators  are  ex¬ 
pected  to  be  substantially  reduced  in  his  jet. 
Accordingly,  no  high  frequency  peaks  are  observed  in  the 
far  field  spectra.  On  the  other  hand-,  the  spectra  meas¬ 
ured  in  the  far  field  of  a  free  jet  by  Lush  (1970)  and 
Ahuja  Alushell  (1973)  exhibit  a  broad  peak.  The  frequency 
of  the  peak  increases  with  angle  to  the  jet  axis,  at  which 
the  far  field  measurements  are  taken.  The  overall  spectra 
are  fairly  broad  due  to  the  random  character  of  the  whole 
turbulent  flow.  It  is  felt  that  with  the  presence  of  the 
plate,  the  flow  is  more  organized  as  seen  from  the  far 
field.  High  frequency  radiators  £ :  .  rirtricted  to  a  short 


distance  near  the  nozzle  by  the  collective  interaction, 
while  downstream,  orderly  structures  dominate  and  radiate 
low  frequency  noise  when  they  impinge  on  the  plate. 

Tam  (1971)  observed  from  shadowgraph  pictures  that 
strong  directional  acoustic  waves  are  emitted  from  the 
shear  layer  close  to  the  exit  of  a  supersonic  free  jet. 
Furthermore,  he  stated  that  "in  a  rather  puzzling  manner 
these  waves  seem  to  exist  only  in  a  limited  region  of 
space  downstream  of  the  nozzle".  Tam  developed  a  theory 
about  the  generation  mechanism  of  these  waves,  based  on 
the  instability  of  the  shear  layer  close  to  the  nozzle  . 
These  observations  and  the  corresponding  theory  are  in 
good  agreement  with  the  two  concepts  presented  in  our 
study,  namely,  the  high  frequency  radiation  by  the  insta¬ 
bility  waves  near  the  nozzle  exit  and  the  collective  in¬ 
teraction  phenomenon. 

In  conclusion,  Fig.  (6-11)  may  help  in  explaining 
the  mechanism  of  far  field  noise  radiation  for  an  imping¬ 
ing  jet.  Based  on  measurements  at  approximately  90  to  the 
jet  axis,  the  far  field  noise  measured  is  a  superposition 
of  low  frequency  and  high  frequency  components.  The  low 
frequency  noise  is  radiated  from  the  plate  by  the  imping¬ 
ing  large  coherent  structures.  The  high  frequency  noise 
radiates  directly  from  the  shear  layer  near  the  nozzle 


CHAPTER  7 


CONCLUDING  REMARKS 

The  present  experimental  investigation  has  success¬ 
fully  Introduced  several  important  concepts  concerning 
both  the  feedback  mechanism  of  a  resonant  impinging  jet 
and  the  far  field  noise  generated  by  jet  impingement  in 
general . 

7.1  THE  FEEDBACK  MECHANISM 

The  feedback  loop  between  the  nozzle  and  the  plate 
consists  of  two  branches:  a  downstream  travelling  wave 
with  a  speed  .62U  due  to  the  convected  large  coherent 
structures,  and  an  upstream  wave  propagating  with  the 
speed  of  sound  of  the  quiescent  medium  in  a  direction  in¬ 
clined  at  an  angle  9  to  the  jet  axis. 

A  phase  lock  necessary  for  self-sustained  oscilla¬ 
tions  is  established  between  the  two  waves  at  the  nozzle 
exit.  Pressure  measurements  (5-1)  and  (5-2))  have 
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indicated  that  the  two  waves  are  in-phase  at  the  nozzle 
exit  plane  during  resonance.  The  phase  lock  led  to  a 
clear  understanding  of  the  resonance  frequency  stages, 
events  which  are  also  observed  in  cavity  and  edge  tones. 
The  number  of  resonant  periods  along  the  entire  feedback 
loop  is  an  Integer  N  and  is  preserved  during  each  stage 
[Fig.  (5-5)].  The  quantum  frequency  jump  between  stages 
takes  place  through  a  quantum  unit  change  in  N. 


It  was  felt  at  the  beginning  that  the  upstream  acous¬ 
tic  waves  would  force  the  initial  shear  layer  instability 
near  the  nozzle.  Instead,  measurements  near  the  nozzle 
exit  have  substantiated  the  idea  that  the  acoustic  forcing 
resulted  in  pulsation  of  the  thin  shear  layer.  This  pul¬ 
sation  forced  a  coherent  merging  of  small  vortices  at  the 
initial  instability  frequency  into  large  vortices  at  the 
resonance  frequency.  The  phenomenon  is  named  "collective 
interaction",  and  its  presence  is  further  supported  by 
similar  observations  in  forced  two-dimensional  free  shear 
layers  by  Ho  and  Huang  (1978). 


A  simple  mathematical  model  was  presented  to  explain 
the  characteristic  behavior  of  measured  correlation  func¬ 
tions  in  a  field  with  two  counter-propagating  waves.  The 
model  provided  evidence  that  the  upstream  wave  is  generat¬ 
ed  by  the  incidence  of  the  c  sr  r earn  travelling  vortices 
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on  the  plate. 

The  analysis  of  the  resonance  phenomenon  in  the  pre- 
aent  study  vas  based  on  near  field  measurements  which  were  ; 
taken  outside  the  jet  column.  The  resulting  picture  in¬ 
cludes  forcing  of  the  shear  layer  by  the  acoustic  waves 
from  the  ambient  side.  The  model  of  Neuworth  (1973)  is 
based  on  acoustic  feedback  from  within  the  jet  core,  but 
unfortunately  no  reliable  measurement  could  be  made  Inside 
a  high  speed  jet.  However,  the  mathematical  analysis  in 
Sec.  5.4  showed  that  the  thin  shear  layer  is  more  suscep¬ 
tible  to  upstream  acoustic  forcing  from  outside  the  jet 
than  from  inside. 


7.2  THE  MECHANISM  OF  NOISE  GENERATION 

The  far  field  noise  for  a  wide  range  of  impinging  jet 

• 

speeds  was  investigated.  Measurements  were  taken  at  90  to 
the  jet  axis  at  the  nozzle  exit  plane.  Surprisingly,  the 
far  field  spectra  not  only  showed  a  low  frequency  peak  but 
also  a  high  frequncy  one  [Fig.  (6-8)]. 

Cross-correlations  between  the  far  field  and  the  near 
field  suggested  that  the  noise  corresponding  to  the  two 
peaks  are  generated  by  two  different  mechanisms,  and  pro¬ 
pagate  along  two  different  paths.  The  low  frequency  noise 

is  mainly  generated  by  the  .rr.ement  of  large  scale 
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structures  on  the  piste  near  r^/d-l.  The  high  frequency 
noise  redletes  directly  from  the  vicinity  of  the  nozzle  i 
exit*  This  noise  is  produced  by  the  high  frequency  lnsta-  j 
bllity  waves  in  the  thin  shear  layer  near  the  nozzle*  j 

It  is  felt  that  the  success  in  localizing  the  above 
regions  of  high  and  low  frequency  noise  radiation  is  due 
to  (i)  the  collective  interaction  which  confines  high  fre¬ 
quency  noise  radiators  to  a  region  near  the  nozzle  exit 
and  (ii)  the  presence  of  the  plate  which  leads  to  feedback 
and  vorticity  stretching,  thereby  enhancing  large  scale 
eddies.  The  impingement  of  these  eddies  on  the  plate 
generates  low  frequency  noise  that  is  radiated  to  the  far 
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2-8)  Block  Diagram  of  Digital  Data  Processing 


(3-l)  Ax i symmetric  Vortical  Rings  in  the  Impinging  Water 
Jet,  Pattern  1:  Small  Spacing  Between  Rings, 
(Re=5000,  x0/d=4) 


Fig.  (3-2)  Pattern  2:  Large  Spacing  Between  Rings,  same 
conditions  as  in  Fig.  (3_l) 


Fig.  (3-5) 


Hydrogen  Bubble  V  i  sua  1  i  z«i  t  ion 
(Re=5000,  x0-4  in.,  d-1.5  in.) 
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(3-8)  Multiple  Exposure  of  Hydrogen  Bubbles  Showing  Shear 

Layer  Growth  and  Ambient  Entrainment,  same  conditions 
as  in  Fig.  (3-7) 
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(ii)  M  =  .9 
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(resonance ) 


Fig.  (3-11) 


Raw  Signals  of  Surface  Pressure  Fluctuations, 

(x0/cM,  ro/d«=0.) 

1 


150 


(3~l4)  Radial  Variation  of  rms  Surface  Pressure 
F?  uctuations 


(3-15)  Variation  of  rms  Surface  Pressure  Fluctuations 


f  rms  Surface  Pressure  Fluctuations 
Locations 


( 3“ 1 7)  Two-Point  Measurements ,  indices 
in  Fig.  (3*19) ,  (M«.8,  xo/d-6) 

(a)  Cross-Spectrum  and  Phase  Angle 

(b)  Coherence  and  Phase  Angle 

(c)  Coherence  and  Phase  Velocity 
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(b)  after  prewhitening 

Fig.  (3-27)  Correlations  between  Surface  Transducers  at  Different 
Azimuthal  Angles,  0Oi"O°  and  002*270°,  (M«,8,  Xq/^, 
roi*ro2"^ • 5d) 


(3-29)  Azimuthal  Cross  Spectrum  (Smooth  Estimate)  and  Phase 
Angle,  same  conditions  as  in  Fig.  (3"27) 
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Fig.  (4-3)  Near  Field  Cross-Spectra  and  Phase  Angle  Functions 
x0/d«5) 


Fig.  (4-5)  Cross-Correlation  of  Near  Field  Pressure  Signals, 
(a)  before  and  (b)  after  Prewhitening,  same  condi ti 
as  in  Fig,  (4-3) 


Radial  Variation  of  rms  Pressure  Fluctuations  in 
Decibels,  at  the  Nozzle  Exit  Plane 
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between  Measured  Cross-Correlations 
the  Model,  (xQ/d*5-5  and  7) 


in  the  Measured  Autocorrelations  and  the  Model 


Fig.  (4-21)  Downstream  and  Upstream  Wave  Propagation  Velocities 
Broadband  Velocity  (0j ,  C?)  and  Phase  Velocity  at 
the  Resonance  Frequency  Cj(fr)  ,  (^(fj.) 
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Measured  data 


'9*  (55)  Resonance  Frequency  Stages, 


a  *  Tam  and  Block  (1978) 
—  Eqn  (5- 4b) 
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rison  between  Cavity  Tones  for  Different  Cavity 
is  ions  as  Measured  by  many  Investigators  (taken 
Tam  and  Block  (1978)  ,  and  Eqn.  (5“4b) 
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( i)  x/d  =.13 ,  r/d  =  .5 


(ii)  x/d  =.92,  r/d  =  1.125 


g.  (5-10a)  Raw  Near  Field  Pressure  Signals  of  a  Non-Resonance 


(ii)  x/d  =  1.31,  r/d  =  I 


Fig.  (5-tOb) 


Raw  Near  Field  Pressure  Signals  of  a  Resonance  Case, 
(M=.9,  x0AM.5 
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Fig.  (6-**)  Cross-Correlations  with  Far  Field  Microphone,  point  2 
Is  at  rQ/d*1 ,  Indices  refer  to  the  sketch  In  Fig.  (*4-2) 
(M-.8,  Xo/d-5.5,  Vd-90) 
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Downstream  convected 
large  scale  structures 


Fig.  (6-6)  Schematic  of  a  Model  for  Far  Field  Noise  Radiation 
due  to  the  Impinging  Large  Scale  Structures 
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Comparison  between  the  Optimum  Time  Delays  of  the 
Correlation  wi th  the  Far  Field  and  the  Model 


between  Far  Field  and  Near  Field  (Close  to 


High  Frequency 
Far -Field  Noise 


Low  Frequency 
Far-Field  Noise 


.  (6-11)  Summary  of  Noise  Generation  Mechanisms  in  an 
Impinging  Jet 
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APPENDIX  A 


STANDING  WAVE  PATTERN 


The  pressure  p  in  Eqn.(4-5)  can  be  written,  without 
1 o 8 s  of  generality,  in  the  form 


a  a  e 


b  e 


(A-») 


The  boundary  condition  at  the  plate  (4-9)  is  satisfied  in 
Eqn,(A-l).  This  equation  can  also  be  written 


far*"  K»  (K~*o)] 

*OC  4bl  *  , 


(A-*) 


where  k«  is  expressed  in  terms  of  its  components  k0  and 

“x 

k2  in  the  x  and  r  directions,  and  the  amplitudes  a  and  b 
r 

are  functions  of  position  .  It  is  assumed  that  the  ampli¬ 


tudes  satisfy  the  inequalities 


u 


C 


Q 


O 


where , 

K  =  “s-  (* - x0 >  +  ( *s  -  r )]  > 

Or\c) 

C  -  r  JL-t  (y  _x0)_  JS£r  (r5-0  J  .  (A-4) 

The  first  tern  in  Eqns.(A-3)  is  a  standing  wave, 
while  the  second  term  represents  a  downstream  propagating 
wave  in  Eqn.(A-3a)  and  an  upstream  propagating  wave  in 
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0 


Eqn . ( A-3  b) •  There  are  nodes  and  antinodes,  but  because  of 
the  second  propagating  terra  the  nodes  are  not  sero  points. 

i 

The  factor  axp  1€  disperses  the  raaxlraun  amplitudes  of  var-  j 

I 

lous  antinode  oscillations  so  they  no  longer  occur  at  the  * 
same  time  (Neuworth  (1973)]. 


Since  observations  indicate  that  resonance  occurs  at 
high  Mach  number  and  small  plate  to  nozzle  distances, £  can 
be  shown  to  be  small.  At  high  M  the  convection  speed  of 
large  scale  structures  is  close  to  the  speed  of  sound. 
Therefore,  at  the  resonance  frequency  the  two  wave  numbers 
kj  and  k2  are  of  the  same  order  of  magnitude ,i .e ., (kj-kj  ) 
is  small.  Also,  (xQ-x)  is  small  at  resonance  and  (rg-r) 
is  generally  small  as  long  as  r  is  in  the  near  field. 


Using  the  expansion 


exp  iC«  1  +  i€-  C2/2  + 
Eqns.(A-3)  simplify  to 
r/d<l .0: 


t'wi 


-  Jibk  e  r  +  (a-b)  t 


i[*ri 


+  o(e) 


r  /d> 1.0; 


Ja  k  e  4.  t  +  0(c) 


which  are  the  forms  used  in  Eqn.  (4-28). 
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